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Astrophysical observations are pointing out huge amounts of "dark matter" and "dark energy" 
needed to explain the observed large scale structure and cosmic dynamics. The emerging picture is 
a spatially flat, homogeneous Universe undergoing the today observed accelerated phase. Despite of 
the good quality of astrophysical surveys, commonly addressed as Precision Cosmology, the nature 
and the nurture of dark energy and dark matter, which should constitute the bulk of cosmological 
- - - matter-energy, are still unknown. Furthermore, up to now, no experimental evidence has been 

found, at fundamental level, to explain such mysterious components. 
, The problem could be completely reversed considering dark matter and dark energy as "short- 

comings" of General Relativity in its simplest formulation (a linear theory in the Ricci scalar R, 
minimally coupled to the standard perfect fluid matter) and claiming for the "correct" theory of 
gravity as that derived by matching the largest number of observational data, without imposing any 
theory a priori. As a working hypothesis, accelerating behavior of cosmic fluid, large scale struc- 
ture, potential of galaxy clusters, rotation curves of spiral galaxies could be reproduced by means of 
extending the standard theory of General Relativity. In other words, gravity could acts in different 
ways at different scales and the above "shortcomings" could be due to incorrect extrapolations of 
the Einstein gravity, actually tested at short scales and low energy regimes. 
' (— I '' After a survey of what is intended for Extended Theories of Gravity in the so called "metric" 

and "Palatini" approaches, we discuss some cosmological and astrophysical applications where the 
I ■ issues related to the dark components are addressed by enlarging the Einstein theory to more 

O ' general f{R) Lagrangians, where f{R) is a generic function of Ricci scalar R, not assumed simply 

linear. Obviously, this is not the flnal answer to the problem of "dark-components" but it can be 
^ . considered as an operative scheme whose aim is to avoid the addition of unknown exotic ingredients 

to the cosmic pie. 

(N ' 

> ■ I. INTRODUCTION 

' General Relativity (GR) is a comprehensive theory of spacetime, gravity and matter. Its formulation implies that 
space and time are not "absolute" entities, as in Classical Mechanics, but dynamical quantities strictly related to the 
distribution of matter and energy. As a consequence, this approach gave rise to a new conception of the Universe 
itself which, for the first time, was considered as a dynamical system. In other words, Cosmology has been enclosed 
' in the realm of Science and not only of Philosophy, as before the Einstein work. On the other hand, the possibility 
, of a scientific investigation of the Universe has led to the formulation of the Standard Cosmological Model [1] which, 
quite nicely, has matched with observations. 

Despite of these results, in the last thirty years, several shortcomings came out in the Einstein theory and people 
^ ■ began to investigate whether GR is the only fundamental theory capable of explaining the gravitational interaction. 
■ Such issues come, essentially, from cosmology and quantum field theory. In the first case, the presence of the Big 
Bang singularity, the flatness and horizon problems [2] led to the statement that Cosmological Standard Model, based 
on the GR and the Standard Model of Particle Physics, is inadequate to describe the Universe at extreme regimes. 
On the other hand, GR is a classical theory which does not work as a fundamental theory, when one wants to achieve 
a full quantum description of spacetime (and then of gravity) . 

Due to these facts and, flrst of all, to the lack of a deflnitive quantum gravity theory, alternative theories have 
been considered in order to attempt, at least, a semi-classical scheme where GR and its positive results could be 
recovered. One of the most fruitful approaches has been that of Extended Theories of Gravity (ETG) which have 
become a sort of paradigm in the study of gravitational interaction. They are based on corrections and enlargements 
of the Einstein theory. The paradigm consists, essentially, in adding higher-order curvature invariants and minimally 
or non- minimally coupled scalar flelds into dynamics which come out from the effective action of quantum gravity [3]. 

Other motivations to modify GR come from the issue of a full recovering of the Mach principle which leads to assume 
a varying gravitational coupling. The principle states that the local inertial frame is determined by some average of 
the motion of distant astronomical objects [4]. This fact implies that the gravitational coupling can be scale-dependent 
and related to some scalar field. As a consequence, the concept of "inertia" and the Equivalence Principle have to be 
revised. For example, the Brans-Dicke theory [5] is a serious attempt to define an alternative theory to the Einstein 
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gravity: it takes into account a variable Newton gravitational coupling, whose dynamics is governed by a scalar field 
non-minimally coupled to the geometry. In such a way, Mach's principle is better implemented [5-7]. 

Besides, every unification scheme as Superstrings, Supergravity or Grand Unified Theories, takes into account 
effective actions where non-minimal couplings to the geometry or higher-order terms in the curvature invariants are 
present. Such contributions are due to one-loop or higher-loop corrections in the high-curvature regimes near the 
full (not yet available) quantum gravity regime [3]. Specifically, this scheme was adopted in order to deal with 
the quantization on curved spacetimes and the result was that the interactions among quantum scalar fields and 
background geometry or the gravitational self-interactions yield corrective terms in the Hilbert-Einstein Lagrangian 
[8]. Moreover, it has been realized that such corrective terms are inescapable in order to obtain the effective action of 
quantum gravity at scales closed to the Planck one [9]. All these approaches are not the "full quantum gravity" but 
are needed as working schemes toward it. 

In summary, higher-order terms in curvature invariants (such as i?,^, R^'^R^y, R'^""^ Rfn^ap , ROR, or RO'^R) or 
non-minimally coupled terms between scalar fields and geometry (such as (p^R) have to be added to the effective 
Lagrangian of gravitational field when quantum corrections are considered. For instance, one can notice that such 
terms occur in the effective Lagrangian of strings or in Kaluza-Klein theories, when the mechanism of dimensional 
reduction is used [10]. 

On the other hand, from a conceptual viewpoint, there are no a priori reason to restrict the gravitational Lagrangian 
to a linear function of the Ricci scalar R, minimally coupled with matter [11]. Furthermore, the idea that there are 
no "exact" laws of physics could be taken into serious account: in such a case, the effective Lagrangians of physical 
interactions are "stochastic" functions. This feature means that the local gauge invariances {i.e. conservation laws) 
are well approximated only in the low energy limit and the fundamental physical constants can vary [12]. 

Besides fundamental physics motivations, all these theories have acquired a huge interest in cosmology due to the 
fact that they "naturally" exhibit inflationary behaviors able to overcome the shortcomings of Cosmological Standard 
Model (based on GR). The related cosmological models seem realistic and capable of matching with the CMBR 
observations [13-15]. Furthermore, it is possible to show that, via conformal transformations, the higher-order and 
non-minimally coupled terms always correspond to the Einstein gravity plus one or more than one minimally coupled 
scalar fields [i6-20]. 

More precisely, higher-order terms appear always as contributions of order two in the field equations. For example, 
a term like R^ gives fourth order equations [21], R OR gives sixth order equations [20, 22], RD'^R gives eighth order 
equations [23] and so on. By a conformal transformation, any 2nd-order derivative term corresponds to a scalar field^: 
for example, fourth-order gravity gives Einstein plus one scalar field, sixth-order gravity gives Einstein plus two scalar 
fields and so on [20, 24]. 

Furthermore, it is possible to show that the /(i?)-gravity is equivalent not only to a scalar-tensor one but also to the 
Einstein theory plus an ideal fiuid [25]. This feature results very interesting if we want to obtain multiple infiationary 
events since an early stage could select "very" large-scale structures (clusters of galaxies today), while a late stage 
could select "small" large-scale structures (galaxies today) [22] . The philosophy is that each inflationary era is related 
to the dynamics of a scalar field. Finafly, these extended schemes could naturafly solve the problem of "graceful exit" 
bypassing the shortcomings of former infiationary models [15, 26]. 

In addition to the revision of Standard Cosmology at early epochs (leading to the Inflation), a new approach is 
necessary also at late epochs. ETGs could play a fundamental role also in this context. In fact, the increasing bulk 
of data that have been accumulated in the last few years have paved the way to the emergence of a new cosmological 
model usually referred to as the Concordance Model. 

The Hubble diagram of Type la Supernovae (hereafter SNela), measured by both the Supernova Cosmology Project 
[27] and the High-z Team [28] up to redshift 2: ~ 1, has been the first evidence that the Universe is undergoing a phase 
of accelerated expansion. On the other hand, balloon born experiments, such as BOOMERanG [29] and MAXIMA 
[30], determined the location of the first and second peak in the anisotropy spectrum of the cosmic microwave 
background radiation (CMBR) strongly pointing out that the geometry of the Universe is spatially flat. If combined 
with constraints coming from galaxy clusters on the matter density parameter I^m, these data indicate that the 
Universe is dominated by a non-clustered fluid with negative pressure, genericafly dubbed dark energy, which is able 
to drive the accelerated expansion. This picture has been further strengthened by the recent precise measurements 
of the CMBR spectrum, due to the WMAP experiment [31-33], and by the extension of the SNela Hubble diagram 
to redshifts higher than 1 [34]. 

After these observational evidences, an overwhelming flood of papers has appeared: they present a great variety 
of models trying to explain this phenomenon. In any case, the simplest explanation is claiming for the well known 
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cosmological constant A [35]. Although it is the best fit to most of the available astrophysical data [31], the ACDM 
model fails in explaining why the inferred value of A is so tiny (120 orders of magnitude lower!) if compared with 
the typical vacuum energy values predicted by particle physics and why its energy density is today comparable to the 
matter density (the so called coincidence problem). 

As a tentative solution, many authors have replaced the cosmological constant with a scalar field rolling down its 
potential and giving rise to the model now referred to as quintessence [36, 37]. Even if successful in fitting the data, 
the quintessence approach to dark energy is still plagued by the coincidence problem since the dark energy and matter 
densities evolve differently and reach comparable values for a very limited portion of the Universe evolution coinciding 
at present era. To be more precise, the quintessence dark energy is tracking matter and evolves in the same way for 
a long time. But then, at late time, somehow it has to change its behavior into no longer tracking the dark matter 
but starting to dominate as a cosmological constant. This is the coincidence problem of quintessence. 

Moreover, it is not clear where this scalar field originates from, thus leaving a great uncertainty on the choice of 
the scalar field potential. The subtle and elusive nature of dark energy has led many authors to look for completely 
different scenarios able to give a quintessential behavior without the need of exotic components. To this aim, it is 
worth stressing that the acceleration of the Universe only claims for a negative pressure dominant component, but 
does not tell anything about the nature and the number of cosmic fiuids filling the Universe. 

This consideration suggests that it could be possible to explain the accelerated expansion by introducing a single 
cosmic fiuid with an equation of state causing it to act like dark matter at high densities and dark energy at low 
densities. An attractive feature of these models, usually referred to as Unified Dark Energy (UDE) or Unified Dark 
Matter (UDM) models, is that such an approach naturally solves, al least plicnomcnologically, the coincidence problem. 
Some interesting examples are the generalized Chaplygin gas [38] , the tachyon field [39] and the condensate cosmology 
[40]. A different class of UDE models has been proposed [41] where a single fiuid is considered: its energy density 
scales with the rcdshift in such a way that the radiation dominated era, the matter era and the accelerating phase 
can be naturally achieved. It is worth noticing that such class of models arc extremely versatile since they can be 
interpreted both in the framework of UDE models and as a two-fiuid scenario with dark matter and scalar field dark 
energy. The main ingredient of the approach is that a generalized equation of state can be always obtained and 
observational data can be fitted. 

Actually, there is still a different way to face the problem of cosmic acceleration. As stressed in [42], it is possible 
that the observed acceleration is not the manifestation of another ingredient in the cosmic pie, but rather the first 
signal of a breakdown of our understanding of the laws of gravitation (in the infra-red hmit). 

From this point of view, it is thus tempting to modify the Fricdmann equations to sec whether it is possible to fit 
the astrophysical data with models comprising only the standard matter. Interesting examples of this kind are the 
Cardassian expansion [43] and the DGP gravity [44]. Moving in this same framework, it is possible to find alternative 
schemes where a quintessential behavior is obtained by taking into account effective models coming from fundamental 
physics giving rise to generalized or higher-order gravity actions [45] (for a comprehensive review see [46]). 

For instance, a cosmological constant term may be recovered as a consequence of a non - vanishing torsion field thus 
leading to a model which is consistent with both SNela Hubble diagram and Sunyaev - Zel'dovich data coming from 
clusters of galaxies [47] . SNela data could also be efficiently fitted including higher-order curvature invariants in the 
gravity Lagrangian [48, 50-52]. It is worth noticing that these alternative models provide naturally a cosmological 
component with negative pressure whose origin is related to the geometry of the Universe thus overcoming the problems 
linked to the physical significance of the scalar field. 

It is evident, from this short overview, the high number of cosmological models which are viable candidates to 
explain the observed accelerated expansion. This abundance of models is, from one hand, the signal of the fact that 
we have a limited number of cosmological tests to discriminate among rival theories, and, from the other hand, that a 
urgent degeneracy problem has to be faced. To this aim, it is useful to remark that both the SNela Hubble diagram 
and the angular size - redshift relation of compact radio sources [53] are distance based methods to probe cosmological 
models so then systematic errors and biases could be iterated. From this point of view, it is interesting to search for 
tests based on time-dependent observables. 

For example, one can take into account the lookback time to distant objects since this quantity can discriminate 
among different cosmological models. The lookback time is observationally estimated as the difference between the 
present day age of the Universe and the age of a given object at redshift z. Such an estimate is possible if the object 
is a galaxy observed in more than one photometric band since its color is determined by its age as a consequence 
of stellar evolution. It is thus possible to get an estimate of the galaxy age by measuring its magnitude in different 
bands and then using stellar evolutionary codes to choose the model that reproduces the observed colors at best. 

Coming to the weak-field-limit approximation, which essentially means considering Solar System scales, ETGs are 
expected to reproduce GR which, in any case, is firmly tested only in this hmit [54]. This fact is matter of debate since 
several relativistic theories do noi reproduce exactly the Einstein results in the Newtonian approximation but, in some 
sense, generalize them. As it was firstly noticed by Stelle [57], a i?^-theory gives rise to Yukawa- like corrections in the 
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Newtonian potential. Such a feature could have interesting physical consequences. For example, some authors claim 
to explain the flat rotation curves of galaxies by using such terms [58] . Others [59] have shown that a conformal theory 
of gravity is nothing else but a fourth-order theory containing such terms in the Newtonian Hmit. Besides, indications 
of an apparent, anomalous, long-range acceleration revealed from the data analysis of Pioneer 10/11, Galileo, and 
Ulysses spacecrafts could be framed in a general theoretical scheme hy taking corrections to the Newtonian potential 
into account [60, 61]. 

In general, any relativistic theory of gravitation yields corrections to the Newton potential (see for example [62]) 
which, in the post-Newtonian (PPN) formalism, could be a test for the same theory [54]. Furthermore the newborn 
gravitational lensing astronomy [63] is giving rise to additional tests of gravity over small, large, and very large scales 
which soon will provide direct measurements for the variation of the Newton coupHng [64], the potential of galaxies, 
clusters of galaxies and several other features of self-gravitating systems. 

Such data will be, very likely, capable of confirming or ruling out the physical consistency of GR or of any ETG. 
In summary, the general features of ETGs are that the Einstein field equations result to be modified in two senses: i) 
geometry can be non- minimally coupled to some scalar field, and/or ii) higher than second order derivative terms in 
the metric come out. In the former case, we generically deal with scalar-tensor theories of gravity; in the latter, we 
deal with higher-order theories. However combinations of non-minimally coupled and higher-order terms can emerge 
as contributions in effective Lagrangians. In this case, we deal with higher-order-scalar-tensor theories of gravity. 

Considering a mathematical viewpoint, the problem of reducing more general theories to Einstein standard form has 
been extensively treated; one can see that, through a "Legendre" transformation on the metric, higher-order theories, 
under suitable regularity conditions on the Lagrangian, take the form of the Einstein one in which a scalar field (or 
more than one) is the source of the gravitational field (see for example [11, 65-67]); on the other side, as discussed 
above, it has been studied the mathematical equivalence between models with variable gravitational coupling with 
the Einstein standard gravity through suitable conformal transformations (sec [68, 69]). 

In any case, the debate on the physical meaning of conformal transformations is far to be solved [see [71] and 
references therein for a comprehensive review]. Several authors claim for a true physical difference between Jordan 
frame (higher-order theories and/or variable gravitational coupHng) since there are experimental and observational 
evidences which point out that the Jordan frame could be suitable to better match solutions with data. Others state 
that the true physical frame is the Einstein one according to the energy theorems [67]. However, the discussion is 
open and no definitive statement has been formulated up to now. 

The problem should be faced from a more general viewpoint and the Palatini approach to gravity could be useful 
to this goal. The Palatini approach in gravitational theories was firstly introduced and analyzed by Einstein himself 
[72]. It was, however, called the Palatini approach as a consequence of an historical misunderstanding [73, 74]. 

The fundamental idea of the Palatini formalism is to consider the (usually torsion-less) connection F, entering the 
definition of the Ricci tensor, to be independent of the metric g defined on the spacetime The Palatini formulation 
for the standard Hilbert-Einstein theory results to be equivalent to the purely metric theory: this follows from the fact 
that the field equations for the connection F, firstly considered to be independent of the metric, give the Levi-Civita 
connection of the metric g. As a consequence, there is no reason to impose the Palatini variational principle in the 
standard Hilbert-Einstein theory instead of the metric variational principle. 

However, the situation completely changes if we consider the ETGs, depending on functions of curvature invariants, 
as f{R), or non-minimally coupled to some scalar field. In these cases, the Palatini and the metric variational principle 
provide different field equations and the theories thus derived differ [67, 75]. The relevance of Palatini approach, in 
this framework, has been recently proven in relation to cosmological applications [45, 46, 76-78]. 

It has also been studied the crucial problem of the Newtonian potential in alternative theories of Gravity and its 
relations with the conformal factor [80]. From a physical viewpoint, considering the metric g and the connection F as 
independent fields means to decouple the metric structure of spacetime and its geodesic structure (being, in general, 
the connection T not the Levi-Civita connection of g). The chronological structure of spacetime is governed by g 
while the trajectories of particles, moving in the spacetime, are governed by F. 

This decoupling enriches the geometric structure of spacetime and generalizes the purely metric formalism. This 
metric-affine structure of spacetime is naturally translated, by means of the same (Palatini) field equations, into a bi- 
metric structure of spacetime. Beside the physical metric g, another metric h is involved. This new metric is related, 
in the case of /(i?)-gravity, to the connection. As a matter of fact, the connection F results to be the Levi-Civita 
connection of h and thus provides the geodesic structure of spacetime. 

If we consider the case of non-minimally coupled interaction in the gravitational Lagrangian (scalar-tensor theories) , 
the new metric h is related to the non-minimal coupling. The new metric h can be thus related to a different geometric 
and physical aspect of the gravitational theory. Thanks to the Palatini formalism, the non-minimal coupling and the 
scalar field, entering the evolution of the gravitational fields, are separated from the metric structure of spacetime. 
The situation mixes when we consider the case of higher-order-scalar-tensor theories. Due to these features, the 
Palatini approach could greatly contribute to clarify the physical meaning of conformal transformation [79]. 
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In this review paper, without claiming for completeness, we want to give a survey on the formal and physical aspects 
of ETGs in metric and Palatini approaches, considering the cosmological and astrophysical applications of some ETG 
models. 

The layout is the following. Sect. II is a rapid overview of GR. We summarize what a good theory of gravity is 
requested to do and what the foundations of the Einstein theory are. The goal is to demonstrate that ETGs have the 
same theoretical bases but, in principle, could avoid some shortcomings of GR which is nothing else but a particular 
case of ETG, f{R) = R. 

The field equations for generic ETGs are derived in Sec. III. Specifically, we discuss two interesting cases: /(-R) and 
scalar-tensor theories considering their relations with GR by conformal transformations. 

The Palatini approach and its intrinsic conformal structure is discussed in Sec. IV giving some peculiar examples. 

Cosmological appHcations are considered in Sec.V. After a short summary of ACDM model, we show that dark 
energy and quintessence issues can be addressed as "curvature effects", if ETGs (in particular /(i?) theories) are 
considered. We work out some cosmological models comparing the solutions with data coming from observational 
surveys. As further result, we show that also the stochastic cosmological background of gravitational waves could be 
"tuned" by ETGs. This fact could open new perspective also in the issues of detection and production of gravitational 
waves which should be investigated not only in the standard framework of GR. 

Sec. VI is devoted to the galactic dynamics under the standard of ETGs. Also in this case, we show that flat rotation 
curves and haloes of spiral galaxies could be explained as curvature effects which give rise to corrections to the Newton 
potential without taking into account huge amounts of dark matter. Discussion and conclusions are drawn in Sec. VII. 

II. WHAT A GOOD THEORY OF GRAVITY HAS TO DO: GENERAL RELATIVITY AND ITS 

EXTENSIONS 

From a phenomenological point of view, there are some minimal requirements that any relativistic theory of gravity 
has to match. First of all, it has to explain the astrophysical observations (e.g. the orbits of planets, the potential of 
self-gravitating structures) . 

This means that it has to reproduce the Newtonian dynamics in the weak-energy limit. Besides, it has to pass the 
classical Solar System tests which are all experimentally well founded [54]. 

As second step, it should reproduce galactic dynamics considering the observed baryonic constituents (e.g. luminous 
components as stars, sub- luminous components as planets, dust and gas), radiation and Newtonian potential which 
is, by assumption, extrapolated to galactic scales. 

Thirdly, it should address the problem of large scale structure (e.g. clustering of galaxies) and finally cosmological 
dynamics, which means to reproduce, in a self-consistent way, the cosmological parameters as the expansion rate, the 
Hubble constant, the density parameter and so on. Observations and experiments, essentially, probe the standard 
baryonic matter, the radiation and an attractive overall interaction, acting at all scales and depending on distance: 
the gravity. 

The simplest theory which try to satisfies the above requirements was formulated by Albert Einstein in the years 
1915-1916 [81] and it is known as the Theory of General Relativity. It is firstly based on the assumption that space 
and time have to be entangled into a single spacetime structure, which, in the limit of no gravitational forces, has to 
reproduce the Minkowski spacetime structure. Einstein profitted also of ideas earlier put forward by Riemann, who 
stated that the Universe should be a curved manifold and that its curvature should be established on the basis of 
astronomical observations [82]. 

In other words, the distribution of matter has to infiuence point by point the local curvature of the spacetime 
structure. The theory, eventually formulated by Einstein in 1915, was strongly based on three assumptions that the 
Physics of Gravitation has to satisfy. 

The "Principle of Relativity" , that amounts to require all frames to be good frames for Physics, so that no preferred 
inertial frame should be chosen a priori (if any exist). 

The "Principle of Equivalence" , that amounts to require inertial effects to be locally indistinguishable from gravi- 
tational effects (in a sense, the equivalence between the inertial and the gravitational mass). 

The "Principle of General Covariance", that requires field equations to be "generally covariant" (today, we would 
better say to be invariant under the action of the group of all spacetime diffeomorphisms) [83] . 

And - on the top of these three principles - the requirement that causality has to be preserved (the " Principle of 
Causality" , i.e. that each point of spacetime should admit a universally valid notion of past, present and future). 

Let us also recall that the older Newtonian theory of spacetime and gravitation - that Einstein wanted to reproduce 
at least in the limit of small gravitational forces (what is called today the "post-Newtonian approximation") - required 
space and time to be absolute entities, particles moving in a preferred inertial frame following curved trajectories, the 
curvature of which (i.e., the acceleration) had to be determined as a function of the sources (i.e., the "forces"). 
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On these bases, Einstein was led to postulate that the gravitational forces have to be expressed by the curvature 
of a metric tensor field ds^ — g^i,dx^dx'' on a four-dimensional spacetime manifold, having the same signature of 
Minkowski metric, i.e., the so-called "Lorentzian signature", herewith assumed to be — , — , — ). He also postulated 
that spacetime is curved in itself and that its curvature is locally determined by the distribution of the sources, i.e. - 
being spacetime a continuum - by the four-dimensional generalization of what in Continuum Mechanics is called the 
"matter stress-energy tensor", i.e. a rank-two (symmetric) tensor T^. 

Once a metric g^^v is given, its curvature is expressed by the Riemann (curvature) tensor 

r/a pa pa i per pa per pa i -t \ 

where the comas are partial derivatives. Its contraction 

ixav — RjjyU 1 (2) 

is the "Ricci tensor" and the scalar 

R = R% = g'^'R^,. (3) 

is called the "scalar curvature" of g^i^. Einstein was led to postulate the following equations for the dynamics of 
gravitational forces 

where k = 8ttG, with c = 1 is a coupling constant. These equations turned out to be physically and mathematically 
unsatisfactory. 

As Hilbert pointed out [83], they were not of a variational origin, i.e. there was no Lagrangian able to reproduce 
them exactly (this is slightly wrong, but this remark is unessential here). Einstein replied that he knew that the 
equations were physically unsatisfactory, since they were contrasting with the continuity equation of any reasonable 
kind of matter. Assuming that matter is given as a perfect fluid, that is 

= (P + P)uiJ.Uu - P9IJ.U (5) 

where u^Ui, is a comoving observer, p is the pressure and p the density of the fluid, then the continuity equation 
requires TJ^ to be covariantly constant, i.e. to satisfy the conservation law 

VT™=0, (6) 

where V denotes the covariant derivative with respect to the metric. 

In fact, it is not true that V^R/^u vanishes (unless i? = 0). Einstein and Hilbert reached independently the 
conclusion that the wrong field equations (4) had to be replaced by the correct ones 

g^. = kt;:1 (7) 

where 

Gf_ii, = R^v — i^diivR (8) 

that is currently called the "Einstein tensor" of g^j^v These equations are both variational and satisfy the conservation 
laws (6) since the following relation holds 

VG^, = 0, (9) 

as a byproduct of the so-called "Bianchi identities" that the curvature tensor of g^,^ has to satisfy [1]. 

The Lagrangian that allows to obtain the field equations (7) is the sum of a "matter Lagrangian" Cm, the variational 
derivative of which is exactly T^, i.e. 



5C^ 

SglJ-v 



7 



and of a "gravitational Lagrangian", currently called the Hilbert-Einstein Lagrangian 

Lhe=9'"'R^.uV^ = RV^, (11) 

where \/—g denotes the square root of the value of the determinant of the metric Qfi^. 

The choice of Hilbert and Einstein was completely arbitrary (as it became clear a few years later), but it was 
certainly the simplest one both from the mathematical and the physical viewpoint. As it was later clarified by Levi- 
Civita in 1919, curvature is not a "purely metric notion" but, rather, a notion related to the "hnear connection" to 
which "parallel transport" and "covariant derivation" refer [84]. 

In a sense, this is the precursor idea of what in the sequel would be called a "gauge theoretical framework" [85], 
after the pioneering work by Cartan in 1925 [86]. But at the time of Einstein, only metric concepts were at hands 
and his solution was the only viable. 

It was later clarified that the three principles of relativity, equivalence and covariance, together with causality, just 
require that the spacctime structure has to be determined by cither one or both of two fields, a Lorentzian metric g 
and a linear connection F, assumed to be torsionless for the sake of simplicity. 

The metric g fixes the causal structure of spacetime (the hght cones) as weh as its metric relations (clocks and 
rods); the connection F fixes the free-fall, i.e. the locally inertial observers. They have, of course, to satisfy a number 
of compatibility relations which amount to require that photons follow null geodesies of F, so that F and g can be 
independent, a priori, but constrained, a posteriori, by some physical restrictions. These, however, do not impose 
that F has necessarily to be the Levi-Civita connection of g [87]. 

This justifies - at least on a purely theoretical basis - the fact that one can envisage the so-called "alternative theories 
of gravitation", that we prefer to call ^'Extended Theories of Gravitation" since their starting points are exactly those 
considered by Einstein and Hilbert: theories in which gravitation is described by either a metric (the so-called "purely 
metric theories"), or by a hnear connection (the so-called "purely afSne theories") or by both fields (the so-cahed 
"mctric-affinc theories", also known as "first order formalism theories"). In these theories, the Lagrangian is a scalar 
density of the curvature invariants constructed out of both g and F. 

The choice (11) is by no means unique and it turns out that the Hilbert-Einstein Lagrangian is in fact the only choice 
that produces an invariant that is linear in second derivatives of the metric (or first derivatives of the connection). 
A Lagrangian that, unfortunately, is rather singular from the Hamiltonian viewpoint, in much than same way as 
Lagrangians, linear in canonical momenta, are rather singular in Classical Mechanics (see e.g. [88]). 

A number of attempts to generahze GR (and unify it to Electromagnetism) along these lines were followed by 
Einstein himself and many others (Eddington, Weyl, Schrodinger, just to quote the main contributors; see, e.g., [89]) 
but they were eventually given up in the fifties of XX Century, mainly because of a number of difficulties related to 
the definitely more complicated structure of a non- linear theory (where by "non-linear" we mean here a theory that 
is based on non- linear invariants of the curvature tensor) , and also because of the new understanding of Physics that 
is currently based on four fundamental forces and requires the more general "gauge framework" to be adopted (see 
[90]). 

Still a number of sporadic investigations about "alternative theories" continued even after 1960 (see [54] and refs. 
quoted therein for a short history) . The search of a coherent quantum theory of gravitation or the belief that gravity 
has to be considered as a sort of low-energy limit of string theories (see, e.g., [91]) - something that we are not wihing 
to enter here in detail - has more or less recently revitalized the idea that there is no reason to follow the simple 
prescription of Einstein and Hilbert and to assume that gravity should be classically governed by a Lagrangian linear 
in the curvature. 

Further curvature invariants or non-linear functions of them should be also considered, especially in view of the 
fact that they have to be included in both the semi-classical expansion of a quantum Lagrangian or in the low-energy 
limit of a string Lagrangian. 

Moreover, it is clear from the recent astrophysical observations and from the current cosmological hypotheses that 
Einstein equations are no longer a good test for gravitation at Solar System, galactic, extra-galactic and cosmic scale, 
unless one does not admit that the matter side of Eqs.(7) contains some kind of exotic matter-energy which is the 
"dark matter" and "dark energy" side of the Universe. 

The idea which we propose here is much simpler. Instead of changing the matter side of Einstein Equations (7) 
in order to fit the "missing matter-energy" content of the currently observed Universe (up to the 95% of the total 
amount!), by adding any sort of inexplicable and strangely behaving matter and energy, we claim that it is simpler and 
more convenient to change the gravitational side of the equations, admitting corrections coming from non-linearities 
in the Lagrangian. However, this is nothing else but a matter of taste and, since it is possible, such an approach 
should be explored. Of course, provided that the Lagrangian can be conveniently tuned up (i.e., chosen in a huge 
family of allowed Lagrangians) on the basis of its best fit with all possible observational tests, at all scales (solar, 
galactic, extragalactic and cosmic). 
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Something that - in spite of some commonly accepted but disguised opinion - can and should be done before rejecting 
a priori a non-linear theory of gravitation (based on a non-singular Lagrangian) and insisting that the Universe has to 
be necessarily described by a rather singular gravitational Lagrangian (one that does not allow a coherent perturbation 
theory from a good Hamiltonian viewpoint) accompanied by matter that does not follow the behavior that standard 
baryonic matter, probed in our laboratories, usually satisfies. 



III. 



THE EXTENDED THEORIES OF GRAVITY 



With the above considerations in mind, let us start with a general class of higher-order-scalar-tensor theories in 
four dimensions ^ given by the action 



A 



-I 



d^xV^ [F{R, DR, D^R, .n'^R, <t>) - -g^" c^;^(t>;u + /:„ 



(12) 



where F is an unspecified function of curvature invariants and of a scalar field (j). The term Cm, as above, is the 
minimally coupled ordinary matter contribution. We shall use physical units SttG = c = = 1; e is a constant which 
specifies the theory. Actually its values can be £ = ±1, fixing the nature and the dynamics of the scalar field which 
can be a standard scalar field, a phantom field or a field without dynamics (see [119, 120] for details). 
In the metric approach, the field equations are obtained by varying (12) with respect to g^v We get 



1 

g 



TM- + ^g'^'^iF - GR) + ig'^^g'"' - g^^'^g^^Q-M 



; dF 
dU^R 



-5" 5 

where G^'^ is the above Einstein tensor and 



dF 
dD'R 



(13) 



(14) 



The differential Eqs.(13) are of order (2fc -|- 4). The stress-energy tensor is due to the kinetic part of the scalar field 
and to the ordinary matter: 



(15) 



The (eventual) contribution of a potential V{(1)) is contained in the definition of F. From now on, we shall indicate by 
a capital F a Lagrangian density containing also the contribution of a potential V{(j)) and by F{(f>), f{R), or f{R, dR) 
a function of such fields without potential. 

By varying with respect to the scalar field </>, we obtain the Klein-Gordon equation 



en(f) ■■ 



dF 



(16) 



Several approaches can be used to deal with such equations. For example, as we said, by a conformal transformation, 
it is possible to reduce an ETG to a (multi) scalar-tensor theory of gravity [18-20, 62, 92]. 
The simplest extension of GR is achieved assuming 



F = f{R), s = 0, 



(17) 



For the aims of this review, we do not need more complicated invariants like Rfi^R^^ , RnvaijR'^'^"^ , C^j^^apC^'^"!^ which are also 
possible. 
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in the action (12); /(i?) is an arbitrary (analytic) function of the Ricci curvature scalar R. We are considering here 
the simplest case of fourth-order gravity but we could construct such kind of theories also using other invariants in 
R^i, or R'^^„- The standard Hilbert-Einstein action is, of course, recovered for f{R) = R. Varying with respect to 
5a/3, we get the field equations 



(18) 



which are fourth-order equations due to the term f'{Ry^"\ the prime indicates the derivative with respect to R. 
Eq.(18) is also the equation for T^i, ~ when the matter term is absent. 
By a suitable manipulation, the above equation can be rewritten as: 



G. 



1 



a/3 



f'{R) {2 



-g^p [f{R) - Rf{R)] + f{R),o,p - g^pUfiR) 



(19) 



where the gravitational contribution due to higher-order terms can be simply reinterpreted as a stress-energy tensor 
contribution. This means that additional and higher-order terms in the gravitational action act, in principle, as a 
stress-energy tensor, related to the form of f{R). Considering also the standard perfect-fiuid matter contribution, we 
have 



1 r 1 1 

= I -5.^ [f{R) - Rf'iR)] + f'{R),^p - 5a/3n/'(J?) j + 



T>curv 
af3 



all 



(20) 



where is an effective stress-energy tensor constructed by the extra curvature terms. In the case of GR, T^'" 

identically vanishes while the standard, minimal coupling is recovered for the matter contribution. The peculiar 
behavior of f{R) = R is due to the particular form of the Lagrangian itself which, even though it is a second order 
Lagrangian, can be non-covariantly rewritten as the sum of a first order Lagrangian plus a pure divergence term. The 
Hilbert-Einstein Lagrangian can be in fact recast as follows: 



Lhe = I^HE\f—g = 



where: 



-gg 



a/3 



dC 

dRat 



(21) 



(22) 



r is the Lcvi-Civita connection of g and il"^;^ is a quantity constructed out with the variation of F [1]. Since i/f^^j is not 
a tensor, the above expression is not covariant; however a standard procedure has been studied to recast covariance 
in the first order theories [93]. This clearly shows that the field equations should consequently be second order and 
the Hilbert-Einstein Lagrangian is thus degenerate. 

From the action (12), it is possible to obtain another interesting case by choosing 



F = F{sP)R-V{cl>), 



-1. 



In this case, we get 



/ 



F{4,)R+-g^''<i>,^4>;.-V{cfy 



(23) 



(24) 



V((p) and F{4)) are generic functions describing respectively the potential and the coupling of a scalar field (j). The 
Brans-Dicke theory of gravity is a particular case of the action (24) for V{(j))=Q [94]. The variation with respect to 
g^i, gives the second-order field equations 



F((/.)G^. = F{cl>) 



1 



Tt.-g^.^^gF{^)+m 



(25) 



here Dg is the d'Alembert operator with respect to the metric g The energy-momentum tensor relative to the scalar 
field is 

T^u = <t>;t^4>;u - \g^..4>;a<i>'^ + g^..V{(t>) (26) 

The variation with respect to (j) provides the Klein - Gordon equation, i.e. the field equation for the scalar field: 

ng(j) - RF^{(1>) + V4,{(t>) = (27) 

where F^ = dF{(f>)/d4), = dV{4))/d([). This last equation is equivalent to the Bianchi contracted identity [95]. 
Standard fiuid matter can be treated as above. 
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A. Conformal transformations 



Let us now introduce conformal transformations to show that any higher-order or scalar-tensor theory, in absence 
of ordinary matter, e.g. a perfect fluid, is conformally equivalent to an Einstein theory plus minimally coupled scalar 
fields. If standard matter is present, conformal transformations allow to transfer non-minimal coupling to the matter 
component [67]. The conformal transformation on the metric g^^ is 

5A.^' = e'^^A.- (28) 
in which e^'^ is the conformal factor. Under this transformation, the Lagrangian in (24) becomes 

(29) 

-6Fa;;„a;.« + ^5""'^;/.'^;. - e-^^V 

in which R and Dg are the Ricci scalar and the d'Alembert operator relative to the metric g. Requiring the theory in 
the metric g^j, to appear as a standard Einstein theory [70], the conformal factor has to be related to F, that is 

g2a. ^ _2p (30) 

where F must be negative in order to restore physical coupling. Using this relation and introducing a new scalar field 
4> and a new potential V , defined respectively by 



2F^ m</')) = ^^, (31) 



we see that the Lagrangian (29) becomes 




-g[FR+ -<?^''</>;^</);. -V]=^g[--R+ -</>;„0;^ - V ) (32) 



which is the usual Hilbert-Einstein Lagrangian plus the standard Lagrangian relative to the scalar field (l>. Therefore, 
every non-minimally coupled scalar-tensor theory, in absence of ordinary matter, e.g. perfect fluid, is conformally 
equivalent to an Einstein theory, being the conformal transformation and the potential suitably deflned by (30) and 
(31). The converse is also true: for a given F(0), such that 3^^^ - F > 0, we can transform a standard Einstein 
theory into a non-minimally coupled scalar-tensor theory. This means that, in principle, if we are able to solve the field 
equations in the framework of the Einstein theory in presence of a scalar field with a given potential, we should be able 
to get the solutions for the scalar-tensor theories, assigned by the coupling F{(t)), via the conformal transformation 
(30) with the constraints given by (31). Following the standard terminology, the "Einstein frame" is the framework of 
the Einstein theory with the minimal coupling and the "Jordan frame" is the framework of the non-minimally coupled 
theory [67]. 

In the context of alternative theories of gravity, as previously discussed, the gravitational contribution to the stress- 
energy tensor of the theory can be reinterpreted by means of a conformal transformation as the stress-energy tensor 
of a suitable scalar field and then as "matter" like terms. Performing the conformal transformation (28) in the field 
equations (19), we get: 

Ga0 = {^5a/3 [f{R) - Rf{R)] + f'{R);a0 ' + (33) 



+2 ( iO-a;l3 + QafdOu! - W;aW;/3 + ^gafjCO-^Ul'^ 



We can then choose the conformal factor to be 



a; = iln|/'(i?)|, (34) 
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which has now to be substituted into (20) . RescaHng w in such a way that 

k(f) = w , 

and k = -^1/6, we obtain the Lagrangian equivalence 



and the Einstein equations in standard form 



V 



with the potential 



2 [P{^)-M{e''"^)e''^] 2 



1 f{R) - Rf'iR) 



f'{Rf 



(35) 



(36) 



(37) 



(38) 



Here N is the inverse function of V^cj)) and 'P((/>) = / exp{2k(f>)dM. However, the problem is completely solved if 
V'id)) can be analytically inverted. In summary, a fourth-order theory is conformally equivalent to the standard 
second-order Einstein theory plus a scalar field (sec also [11, 65]). 

This procedure can be extended to more general theories. If the theory is assumed to be higher than fourth order, we 
may have Lagrangian densities of the form [20, 73], 

C = jC{R,nR,..n''R). (39) 

Every □ operator introduces two further terms of derivation into the field equations. For example a theory like 

C = RnR, (40) 

is a sixth-order theory and the above approach can be pursued by considering a conformal factor of the form 

1 



c.= -ln 



dC ^ dC 
+ 0- 



dR dUR 



(41) 



In general, increasing two orders of derivation in the field equations {i.e. for every term Di?), corresponds to adding 
a scalar field in the conformally transformed frame [20]. A sixth-order theory can be reduced to an Einstein theory 
with two minimally coupled scalar fields; a 2n-order theory can be, in principle, reduced to an Einstein theory plus 
(n— l)-scalar fields. On the other hand, these considerations can be directly generalized to higher-order-scalar-tensor 
theories in any number of dimensions as shown in [17]. 

As concluding remarks, we can say that conformal transformations work at three levels: i) on the Lagrangian 
of the given theory; ii) on the field equations; in) on the solutions. The table below summarizes the situation for 
fourth-order gravity (FOG), non-minimally coupled scalar-tensor theories (NMC) and standard Hilbort-Einstein (HE) 
theory. Clearly, direct and inverse transformations correlate all the steps of the table but no absolute criterion, at 
this point of the discussion, is able to select which is the "physical" framework since, at least from a mathematical 
point of view, all the frames are equivalent [67]. This point is up to now unsolved even if wide discussions are present 
in literature [71]. 



C-FOG 

I 

FOG Eqs. 

I 

FOG Solutions 



^NMC 
I 

NMC Eqs. 

I 

NMC Solutions 



I 

Einstein Eqs. 

I 

Einstein Solutions 



IV. THE PALATINI APPROACH AND THE INTRINSIC CONFORMAL STRUCTURE 



As we said, the Palatini approach, considering g and F as independent fields, is "intrinsically" bi-metric and capable 
of disentangling the geodesic structure from the chronological structure of a given manifold. Starting from these 
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considerations, conformal transformations assume a fundamental role in defining the affine connection which is merely 
"Levi-Civita" only for the Hilbert-Einstein theory. 

In this section, we work out examples showing how conformal transformations assume a fundamental physical role 
in relation to the Palatini approach in ETGs [79]. 

Let us start from the case of fourth-order gravity where Palatini variational principle is straightforward in showing 
the differences with Hilbert-Einstein variational principle, involving only metric. Besides, cosmological applications 
of f{R) gravity have shown the relevance of Palatini formalism, giving physically interesting results with singularity 
- free solutions [76]. This last nice feature is not present in the standard metric approach. 

An important remark is in order at this point. The Ricci scalar entering in ,f{B) is B = R{g,T) = g"^ B.ap{T) 
that is a generalized Ricci scalar and i?^i^(r) is the Ricci tensor of a torsion-less connection F, which, a priori, has 
no relations with the metric g of spacetime. The gravitational part of the Lagrangian is controlled by a given real 
analytical function of one real variable f{R), while ^—g denotes a related scalar density of weight 1. Field equations, 
deriving from the Palatini variational principle are: 

/'(ii)%.)(F) - \f{R)g^^ = (42) 

'^l{V^f'{R)9n = (43) 

where V'" is the covariant derivative with respect to F. As above, we assume 8ttG = 1. We shall use the standard 
notation denoting by R(ij,u) the symmetric part of R^v, i-e. R{tj,v) = \{Rijiv + Run)- 

In order to get (43), one has to additionally assume that is functionally independent of F; however it may contain 

g 

metric covariant derivatives V of fields. This means that the matter stress-energy tensor = TJ^{g, ^) depends 
on the metric g and some matter fields denoted here by 5", together with their derivatives (covariant derivatives with 
respect to the Levi-Civita connection of g). From (43) one sees that \/—gf'{R)g'^^ is a symmetric twice contravariant 
tensor density of weight 1. As previously discussed in [75, 79], this naturally leads to define a new metric hf^v, such 
that the following relation holds: 

V^nR)g^- = ^h'^" . (44) 

This ansatz is suitably made in order to impose F to be the Levi-Civita connection of h and the only restriction is 
that \/—gf'{R)g'^^ should be non-degenerate. In the case of Hilbert-Einstein Lagrangian, it is f'{R) = 1 and the 
statement is trivial. 

The above Eq.(44) imposes that the two metrics h and g are conformally equivalent. The corresponding conformal 
factor can be easily found to be f'{R) (in dim Ai — -i) and the conformal transformation results to be ruled by: 

= f'{R)9i.. (45) 

Therefore, as it is well known, Eq.(43) implies that F = TLc(h) and i?(^^)(r) = R^v{h) = Rf^u- Field equations can 
be supplemented by the scalar- valued equation obtained by taking the trace of (42) , (we define r = trT) 

fiR)R - 2f{R) = g^f'T^^ = (46) 

which controls solutions of (43). We shall refer to this scalar- valued equation as the structural equation of the 
spacetime. In the vacuum case (and spacetimes filled with radiation, such that r™ = 0) this scalar-valued equation 
admits constant solutions, which are different from zero only if one add a cosmological constant. This means that the 
universality of Einstein field equations holds [75], corresponding to a theory with cosmological constant [96]. 

In the case of interaction with matter fields, the structural equation (45), if explicitly solvable, provides an expression 
of i? = F{t), where F is a generic function, and consequently both f{R) and f'{R) can be expressed in terms of r. 
The matter content of spacetime thus rules the bi-metric structure of spacetime and, consequently, both the geodesic 
and metric structures which are intrinsically different. This behavior generalizes the vacuum case and corresponds to 
the case of a time-varying cosmological constant. In other words, due to these features, conformal transformations, 
which allow to pass from a metric structure to another one, acquire an intrinsic physical meaning since "select" metric 
and geodesic structures which, for a given ETG, in principle, do not coincide. 

Let us now try to extend the above formalism to the case of non-minimally coupled scalar-tensor theories. The 
effort is to understand if and how the bi-metric structure of spacetime behaves in this cases and which could be its 
geometric and physical interpretation. 

We start by considering scalar-tensor theories in the Palatini formalism, calling Ai the action functional. After, we 
take into account the case of decoupled non-minimal interaction between a scalar-tensor theory and a f{R) theory. 
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calling A2 this action functional. We finally consider the case of non-minimal-couplcd interaction between the scalar 
field (p and the gravitational fields ((7,r), caUing A3 the corresponding action functional. Particularly significant is, 
in this case, the limit of low curvature R. This resembles the physical relevant case of present values of curvatures of 
the Universe and it is important for cosmological applications. 

The action (24) for scalar-tensor gravity can be generalized, in order to better develop the Palatini approach, as: 

Ai= [F{<P)R +i^^cj>\7%- V{cP) + £„(*, V ^d^x . (47) 

As above, the values of e = ±1 selects between standard scalar field theories and quintessence (phantom) field theories. 
The relative "signature" can be selected by conformal transformations. Field equations for the gravitational part of 
the action are, respectively for the metric g and the connection F: 



r F(0)[i?(^,) - ^Rg^,] = T^t + 
\'7li^F{^)g^-) = 

is the same defined in (42). For matter fields we have the following field equations: 

eU(t) = -V^{(t)) + F^{<j>)R 



SCr, 



= 



(48) 



(49) 



(5* 

In this case, the structural equation of spacetime implies that: 

which expresses the value of the Ricci scalar curvature in terms of the traces of the stress-energy tensors of standard 
matter and scalar field (we have to require -^(0) 7^ 0). The bi-metric structure of spacetime is thus defined by the 
ansatz: 

^/=^F(.^)5''^ = sT^h''^ (51) 
such that g and h result to be conformally related 

V = F{^)g^, . (52) 

The conformal factor is exactly the interaction factor. From (50), it follows that in the vacuum case r"^ = and 
j.m _ Q. -(-j-^jg theorj- is equivalent to the standard Einstein one without matter. On the other hand, for F{(f>) = Fq we 
recover the Einstein theory plus a minimally coupled scalar field: this means that the Palatini approach intrinsically 
gives rise to the conformal structure (52) of the theory which is trivial in the Einstein, minimally coupled case. 

As a further step, let us generaHze the previous results considering the case of a non-minimal coupling in the 
framework of f{R) theories. The action functional can be written as: 



/ 



-g [F{<p)f{R) + is/^cj>V <P- V{cP) + £„(*, V ^)]d''x (53) 



where f{R) is, as usual, any analytical function of the Ricci scalar R. Field equations (in the Palatini formalism) for 
the gravitational part of the action are: 

r f (0)[/'(i?)i?(^.) - i/(i?)5M.] = T^L + ,.4^ 
\vl{V^F{^)f'{R)gn = 0- ^ ' 

For scalar and matter fields we have, otherwise, the following field equations: 

(enct>=-V^{cl>) + ^F4ct>)f{R) 
i Mm _ n 

where the non-minimal interaction term enters into the modified Klein-Gordon equations. In this case the structural 
equation of spacetime impHes that: 

_0 _|_ _m 

f{R)R-2f{R) = ^^. (56) 
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We remark again that this equation, if solved, expresses the value of the Ricci scalar curvature in terms of traces of 
the stress-energy tensors of standard matter and scalar field (we have to require again that -F(0) ^ 0). The bi- metric 
structure of spacetime is thus defined by the ansatz: 



such that g and h result to be conformally related by: 

K. = F{<l>)f{R)g^,. 



(57) 



(58) 



Once the structural equation is solved, the conformal factor depends on the values of the matter fields {<h, '5) or, more 
precisely, on the traces of the stress-energy tensors and the value of </>. From equation (56), it follows that in the 
vacuum case, i.e. both r"^ = and r™ = 0, the universality of Einstein field equations still holds as in the case of 
minimally interacting f{R) theories [75]. The validity of this property is related to the decoupHng of the scalar field 
and the gravitational field. 

Let us finally consider the case where the gravitational Lagrangian is a general function of (f) and R. The action 
functional can thus be written as: 



A. = 



y(</.)+£„(*,V*)]d*a; 



Field equations for the gravitational part of the action are: 



I K{4>,R) 



OR 



(^ 



R 



OKU 



M,R)g, 



dR 



0. 



rp(^ _|_ rpm 

1X1/ * ixv 



(59) 



(60) 



For matter fields, we have: 




The structural equation of spacetime can be expressed as: 



dK{(i), R) 
dR 



R - 2K{(j}, R) = t'^+ r" 



(61) 



(62) 



This equation, if solved, expresses again the form of the Ricci scalar curvature in terms of traces of the stress-energy 
tensors of matter and scalar field (we have to impose regularity conditions and, for example, K{(f),R) ^ 0). The 
bi-metric structure of spacetime is thus defined by the ansatz: 



-dK{(t>, R) 



dR 



-hh'"' 



such that g and h result to be conformally related by 



h 



dK{(j), R ) 
OR 



(63) 



(64) 



Again, once the structural equation is solved, the conformal factor depends just on the values of the matter fields and 
(the trace of) their stress energy tensors. In other words, the evolution, the definition of the conformal factor and the 
bi-metric structure is ruled by the values of traces of the stress-energy tensors and by the value of the scalar field (j). 
In this case, the universality of Einstein field equations does not hold anymore in general. This is evident from (62) 
where the strong coupHng between R and (j) avoids the possibility, also in the vacuum case, to achieve simple constant 
solutions. 

We consider, furthermore, the case of small values of R, corresponding to small curvature spacetimes. This limit 
represents, as a good approximation, the present epoch of the observed Universe under suitably regularity conditions. 
A Taylor expansion of the analytical function K[(j), R) can be performed: 



K{cj>, R) = Ko{cf>) + Ki{cj>)R + o{R^) 



(65) 
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where only the first leading term in R is considered and we have defined: 

( Ko{cP)=K{(^,R)r^o 

Substituting this expression in (62) and (64) we get (neglecting higher order approximations in R) the structural 
equation and the bi-metric structure in this particular case. From the structural equation, we get: 

such that the value of the Ricci scalar is always determined, in this first order approximation, in terms of t*, t™, (p. 
The bi-metric structure is, otherwise, simply defined by means of the first term of the Taylor expansion, which is 

= Ki{(l))g^^ . (68) 

It reproduces, as expected, the scalar-tensor case (52). In other words, scalar-tensor theories can be recovered in a 
first order approximation of a general theory where gravity and non-minimal couplings are any (compare (67) with 
(56)). This fact agrees with the above considerations where Lagrangians of physical interactions can be considered as 
stochastic functions with local gauge invariance properties [12]. 

Finally we have to say that there are also bi-metric theories which cannot be conformally related (see for example 
the summary of alternative theories given in [54]) and torsion field should be taken into account, if one wants to 
consider the most general viewpoint [55, 56]. We will not take into account these general theories in this review. 

After this short review of ETGs in metric and Palatini approach, we are going to face some remarkable applications 
to cosmology and astrophysics. In particular, we deal with the straightforward generalization of GR, the f{R) gravity, 
showing that, in principle, no further ingredient, a part a generalized gravity, could be necessary to address issues as 
missing matter (dark matter) and cosmic acceleration (dark energy) . However what we are going to consider here are 
nothing else but toy models which are not able to fit the whole expansion history, the structure growth law and the 
CMB anisotropy and polarization. These issues require more detailed theories which, up to now, are not available 
but what we are discussing could be a useful working paradigm as soon as refined experimental tests to probe such 
theories will be proposed and pursued. In particular, we wih outline an independent test, based on the stochastic 
background of gravitational waves, which could be extremely useful to discriminate between ETGs and GR or among 
the ETGs themselves. In this latter case, the data delivered from ground-based interferometers, like VIRGO and 
LIGO, or the forthcoming space interferometer LISA, could be of extreme relevance in such a discrimination. 

Finally, we do not take into account the well known inflationary models based on ETGs (e.g. [13, 15]) since we 
want to show that also the last cosmological epochs, directly related to the so called Precision Cosmology, can be 
framed in such a new "economic" scheme. 



V. APPLICATIONS TO COSMOLOGY 



As discussed in the Introduction, many rival theories have been advocated to fit the observed accelerated expansion 
and to solve the puzzle of dark energy. 

As a simple classification scheme, we may divide the different cosmological models in three wide classes. According 
to the models of the first class, the dark energy is a new ingredient of the cosmic Hubble fiow, the simplest case being 
the ACDM scenario and its quintessential generalization (the QCDM models). 

This is in sharp contrast with the assumption of UDE models (the second class) where there is a single fluid 
described by an equation of state comprehensive of all regimes of cosmic evolution [41] (the parametric density models 
or generalized Eoff^ models). 

Finally, according to the third class of models, accelerated expansion is the first evidence of a breakdown of the 
Einstein GR (and thus of the Friedmann equations) which has to be considered as a particular case of a more general 
theory of gravity. As an example of this kind of models, we wih consider the /(i?) - gravity [45, 46, 48, 50]. 

Far from being exhaustive, considering these three classes of models allow to explore very different scenarios proposed 
to explain the observed cosmic acceleration [97, 98, 105]. However, from the above considerations, it is possible to 



^ EoS for Equation of State. 
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show that one of the simplest extension of GR, the /(i?) gravity can, in principle, address the dark energy issues both 
in metric and Palatini approach. In this section, without claiming for completeness, we sketch some f{R) models 
matching solutions against some sets of data. The goal is to show that the dark energy issue could be addressed as a 
curvature effect in ETGs. 



A. The ACDM model: the paradigm 

Cosmological constant A has become a textbook candidate to drive the accelerated expansion of the spatially flat 
Universe. Despite its conceptual problems, the ACDM model turns out to be the best fit to a combined analysis 
of completely different astrophysical data ranging from SNela to CMBR anisotropy spectrum and galaxy clustering 
[31, 106]. As a simple generalization, one may consider the QCDM scenario in which the barotropic factor w = p/p 
takes at a certain epoch a negative value with w = —1 corresponding to the standard cosmological constant. Testing 
whether such a barotropic factor deviate or not from —1 is one of the main issue of modern observational cosmology. 
How such a negative pressure fluid drives the cosmic acceleration may be easily understood by looking at the Friedmann 
equations : 

H'^i^^y = ^{Pm+PA) , (69) 

2^+H^ = -PA = -wpA , (70) 

where a{t) is the scale factor of the Universe, the dot denotes differentiation with respect to cosmic time t, H is the 
Hubble parameter and the Universe is assumed spatially flat as suggested by the position of the flrst peak in the 
CMBR anisotropy spectrum (see also Fig.l). )[29-31]. 

From the continuity equation, p + 3i/(/9 + p) — Q, we get for the i-th fluid with pi = WiPi : 

Qi = Oi,oa-'^'+"*^ = Oi,o(l + , (71) 

where 2 = l/a— lis the redshift, Hi = pi/pcrit is the density parameter for the i-th fluid in terms of the critical 
density which, defined in standard units, is pcriL = SHq/SttG and, hereafter, we label all the quantities evaluated 
today with a subscript 0. It is important to stress that Eq.(71) works only for Wi = constant. Inserting this result 
into Eq.(69), one gets : 

H{z) = Ho^^MfiCi- + + ^^A,o(l + ^)3(i+«') . (72) 

The subscript M means all the matter content, inclusive of dark and baryonic components. Using Eqs.(69), (70) and 
the definition of the deceleration parameter q = —aa/a^, one finds : 

1 3 

10= 2 + 2^^^ ~ ■ ^"^^^ 

The SNela Hubble diagram, the large scale galaxy clustering and the CMBR anisotropy spectrum can all be fitted 
by the ACDM model with (flMfii^A) — (0.3,0.7) thus giving go — —0.55, i.e. the Universe turns out to be in an 
accelerated expansion phase. The simplicity of the model and its capability of fitting the most of the data are the 
reasons why the ACDM scenario is the leading candidate to explain the dark energy cosmology. Nonetheless, its 
generalization, QCDM models, i.e. mechanisms allowing the evolution of A from the past are invoked to remove the 
A-problem and the coincidence problem. 

Here, we want to show that assuming f{R) gravity, not strictly linear in R as GR, it is possible to give rise to 
the evolution of the barotropic factor w = p/p, today leading to the value w = —1, and to obtain models capable of 
matching with the observations. However, also if the paradigm could result valid, it is very difficult to address, in the 
same comprehensive f{R) model, different issues as structure formation, nucleosynthesis, Hubble diagram, radiation 
and matter dominated behaviors as we shall discuss below. 

Before considering specific f{R) theories, let us discuss methods to constrain models by samples of data. 
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Figure 1: The CMBR anisotropy spectrum for different values of w. Data points are the WMAP measurements and the best 
fit is obtained for w ~ —1. If w 7^ — 1 the clustering of dark energy has been considered in this plot. 



B. Methods to constrain models by distance and time indicators 



In principle, cosmological models can be constrained using suitable distance and/or time indicators. As a general 
remark, solutions coming from cosmological models have to be matched with observations by using the redshift z as 
the natural time variable for the Hubble parameter, i.e. 



Data can be obtained for various values of redshift z: for example, CMB probes recombination at z ~ 1100 and 
2; ~ 1 via the late integrated Sachs- Wolfe effect; for 10 < z < 100 the planned 21cm observations could give detailed 
information [134]; futuristic LSS surveys and SNe could probe the Universe up to z ~ 4. The method consists in 
building up a reasonable patchwork of data coming from different epochs and then matching them with the same 
cosmological solution ranging, in principle, from inflation to present accelerated era. 

In order to constrain the parameters characterizing the cosmological solution, a reasonable approach is to maximize 
the following likelihood function : 



C cx exp 



xHpY 



(75) 



where p are the parameters characterizing the cosmological solution. The merit function can be deflned as : 
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x^(p) = E 



i=l 



obs 



n{p) - 1.716 



0.062 



^(p) - 0.469 



0.017 



(76) 



Terms entering Eq.(76) can be characterized as follows. For example, the dimensionless coordinate distances y to 
objects at redshifts z are considered in the first term. They are defined as : 



y{z) 



/" 

Jo 



dz' 



(77) 



where E{z) = H{z)/Hq is the normalized Hubble parameter. This is the main quantity which allows to compare the 
theoretical results with data. The function y is related to the luminosity distance Dl = (1 + z)y{z). 

A sample of data at y{z) for the 157 SNela is discussed in the Riess et al. [34] Gold dataset and 20 radio-galaxies are 
in [108]. These authors fit with good accuracy the linear Hubble law at low redshift {z < 0.1) obtaining the Hubble 
dimensionless parameter h = 0.664±0.008 , with h the Hubble constant in units of 100 km Mpc~^. Such a number 
can be consistently taken into account at low redshift. This value is in agreement with Hq = 72±8 km Mpc""*^ 
given by the HST Key project [109] based on the local distance ladder and estimates coming from time delays in 
multiply imaged quasars [110] and Sunyaev - Zel'dovich effect in X-ray emitting clusters [111]. The second term in 
Eq.(76) allows to extend the 2-range to probe y{z) up to the last scattering surface {z > 1000). The shift parameter 
[112, 113] TZ = \/flMy{zis) can be determined from the CMBR anisotropy spectrum, where zis is the rodsliift of the 
last scattering surface which can be approximated as z^ = 1048 (l + 0.00124^^"" '^^^) (1 + ffiw^) with LOi = VLiK^ 
(with i = b,M for baryons and total matter respectively) and {gi,g2) given in [114]. The parameter uib is constrained 
by the baryogenesis calculations contrasted to the observed abundances of primordial elements. Using this method, 
the value ojb = 0.0214±0.0020 is found [115]. 

In any case, it is worth noticing that the exact value of zig has a negligible impact on the results and setting 
zis = 1100 does not change constraints and priors on the other parameters of the given model. The third term in the 
function takes into account the acoustic peak of the large scale correlation function at 100 Mpc separation, 
detected by using 46748 luminous red galaxies (LRG) selected from the SDSS Main Sample [116, 117]. The quantity 



A 



Zlrg 



Zlrg 
E{zlrg) 



y'^iZLRc) 



1/3 



(78) 



is related to the position of acoustic peak where Zlrg = 0.35 is the effective redshift of the above sample. The 
parameter A depends on the dimensionless coordinate distance (and thus on the integrated expansion rate), on VIm 
and E{z). This dependence removes some of the degeneracies intrinsic in distance fitting methods. 

Due to this reason, it is particularly interesting to include ^ as a further constraint on the model parameters using 
its measured value A = 0.469±0.017 [116]. Note that, although similar to the usual introduced in statistics, the 
reduced (i-e., the ratio between the x^ and the number of degrees of freedom) is not forced to be 1 for the best fit 
model because of the presence of the priors on TZ and A and since the uncertainties CTj are not Gaussian distributed, 
but take care of both statistical errors and systematic uncertainties. With the definition (75) of the likelihood function, 
the best fit model parameters are those that maximize C{p). 

In order to implement the above sketched method, much attention, has been devoted to standard candles, i.e. astro- 
physical objects whose absolute magnitude M is known (or may be exactly predicted) a priori so that a measurement 
of the apparent magnitude m immediately gives the distance modulus jj, = m — M . Specifically, the distance to the 
object, estimated in Mpc, is : 



IJ.{z) = 5logDL{z)/Mpc + 25 



(79) 



with Dl{z) the luminosity distance and z the redshift of the object. The number 25 depends on the distance modulus 
calculated in Mpc. The relation between /i and z is what is referred to as Hubble diagram and it is an open window 
on the cosmography of the Universe. Furthermore, the Hubble diagram is a powerful cosmological test since the 
luminosity distance is determined by the expansion rate as : 



Dl{z)^^{1 + z) f 



dz' 

W) 



(80) 
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with E{z) defined above. Being the Hubble diagram related to the luminosity distance and being Dl determined by 
the expansion rate H{z), it is clear why it may be used as an efficient tool to test cosmological models and constrain 
their parameters. 

To this aim, however, it is mandatory that the relation ^ = fj,{z) is measured up to high enough redshift since, 
for low z, Dl reduces to a linear function of the redshift (thus recovering the Hubble law) whatever the background 
cosmological model is. This necessity claims for standard candles that are bright enough to be visible at such high 
redshift that the Hubble diagram may discriminate among different rival theories. SNela are, up to now, the objects 
that best match these requirements. 

It is thus not surprising that the first evidences of an accelerating Universe came from the SNela Hubble diagram 
and dedicated survey (like the SNAP satellite [118]) have been planned in order to increase the number of SNela 
observed and the redshift range probed. 

A reliable compilation of SNela is the Gold dataset released by Ricss et al. [34]. The authors have compiled a 
catalog containing 157 SNela with z in the range (0.01, 1.70) and visual absorption Ay < 0.5. The distance modulus 
of each object has been evaluated by using a set of calibrated methods so that the sample is homogenous in the sense 
that all the SNela have been re-analyzed using the same technique in such a way that the resulting Hubble diagram 
is indeed reliable and accurate. Given a cosmological model assigned by a set of parameters p = (pi, . . . ,Pn), the 
luminosity distance may be evaluated with Eq.(80) and the predicted Hubble diagram contrasted with the observed 
SNela one. Constraints on the model parameters may then be extracted by mean of a - based analysis defining, in 
this case, the above as : 



where (Tj is the error on the distance modulus at redshift Zi and the sum is over the NsNeia SNela observed. It is 
worth stressing that the uncertainty on measurements also takes into account errors on the redshifts and they are not 
Gaussian distributed. 

As a consequence, the reduced (i-e., XsNeia divided by the number of degrees of freedom) for the best fit model 
is not forced to be close to unity. Nonetheless, different models may still be compared on the basis of the value : 
the lower is X^Neia^ the better the model fits the SNela Hubble diagram. 

The method outlined is a simple and quite efficient way to test whether a given model is a viable candidate to 
describe the late time evolution of the Universe. Nonetheless, it is affected by some degeneracies that could be only 
partially broken by increasing the sample size and extending the probed redshift range. A straightforward example may 
help in elucidating this point. Let us consider the fiat concordance cosmological model with matter and cosmological 
constant. It is : 



so that XsNeia '^iH oiily depend on the Hubble constant Hq and the matter density parameter ^Im- Actually, we could 
split the matter term in a baryonic and a non-baryonic part denoting with fli, the baryon density parameter. Since 
both baryons and non baryonic dark matter scales as (1 + ^;)^, E{z) and thus the luminosity distance will depend only 
on the total matter density parameter and we could never constrain fib by fitting the SNela Hubble diagram. Similar 
degeneracies could also happen with other cosmological models thus stressing the need for complementary probes to 
be combined with the SNela data. For a review, sec the contribution by Bob Nichols in this volume. 

To this aim, a recently proposed test, based on the gas mass fraction in galaxy clusters, can be considered. We 
briefiy outHne here the method referring the interested reader to the literature for further details [121, 122]. Both 
theoretical arguments and numerical simulations predict that the baryonic mass fraction in the largest relaxed galaxy 
clusters should be invariant with the redshift (sec, e.g., Rcf. [124]). 

However, this will only appear to be the case when the reference cosmology in making the baryonic mass fraction 
measurements matches the true underlying cosmology. From the observational point of view, it is worth noticing 
that the baryonic content in galaxy clusters is dominated by the hot X - ray emitting intra-cluster gas so that what is 
actually measured is the gas mass fraction fgas and it is this quantity that should be invariant with the redshift within 
the caveat quoted above. Moreover, it is expected that the baryonic fraction in clusters equals the universal ratio 
flb/flM so that fgas should indeed be given by bxflb/^M where the multipHcative factor b is motivated by simulations 
that suggest that the gas fraction is slightly lower than the universal ratio because of processes that convert part of 
the gas into stars or eject it outside the cluster itself. 

Following Ref. [123], we adopt the SCDM model (i.e., a fiat Universe with CIm = 1 and h = 0.5, being h the Hubble 
constant in units of 100 km s^^ Mpc^^) as reference cosmology in making the measurements so that the theoretical 
expectation for the apparent variation of fgas with the redshift is [123] : 




(81) 



E^z) = nM{i + zf + {i-nM) 
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bfth 



M 



nSCDM 



1.5 



(82) 



where D^^'^ and D^°'^ is the angular diameter distance for the SCDM and the model to be tested respectively. 
Da{z) may be evaluated from the luminosity distance Dl{z) as : 



DA{z) = {l + zr^DL{z) (83) 

with Dl{z) given by Eq.(80) above. 

In[123], it has been extensively analyzed the set of simulations in Ref. [124] to get b = 0.824±0.089. For values 
in the la range quoted above, the main results are independent on b. It is worth noticing that, while the angular 
diameter distance depends on E{z) and thus on h and r^M, the prefactor in Eq.(82) makes fgas explicitly depending 
on flb/^M so that a direct estimate of fib is (in principle) possible. Actually, for the models which we are going to 
consider, the quantity that is constrained by the data is the ratio Clh/flM rather than O5 itself. 

To simultaneously take into account both the fit to the SNela Hubble diagram and the test on the fgas data, it is 
convenient to perform a likelihood analysis defining the following likelihood function : 



£(p) oc exp 



X'(P)" 



(84) 



with : 



X =XSNeIa+Xgas+[ p pg 1 +( ^ ) (85) 



where it is possible to define : 



Xgas 



(86) 



i=l 

being fgaa{zi) the measured gas fraction in a galaxy clusters at redshift Zi with an error agi and the sum is over the 
Ngas clusters considered. In order to avoid possible systematic errors in the fgas measurement, it is desirable that 
the cluster is both highly luminous (so that the S/N ratio is high) and relaxed so that both merging processes and 
cooling fiows arc absent. A catalog of 26 large relaxed clusters, with a measurement of both the gas mass fraction 
fgas and the redshift z is in [123]. These data can be used to perform a suitable likelihood analysis. 

Note that, in Eq.(85), we have explicitly introduced two Gaussian priors to better constrain the model parameters. 
First, there is a prior on the Hubble constant h determined by the results of the HST Key project [125] from an 
accurate calibration of a set of different local distance estimators. Further, we impose a constraint on the ratio 
flb/flM by considering the estimates of flbh^ and flMh^ obtained by Tegmark et al. [126] from a combined fit to 
the SNela Hubble diagram, the CMBR anisotropy spectrum measured by WMAP and the matter power spectrum 
extracted from over 200000 galaxies observed by the SDSS collaboration. It is worth noticing that, while our prior 
on h is the same as that used by many authors when applying the fgas test [122, 123], it is common to put a second 
prior on Ctb rather than Ctb/f^M- Actually, this choice can be motivated by the peculiar features of the models which 
one is going to consider. 

With the definition (84) of the likelihood function, the best fit model parameters are those that maximize C{p). 
However, to constrain a given parameter Pi, one resorts to the marginalized likelihood function defined as : 



^PiiPi) oc j dpi... j dpi-i j dpi+i ... j dpnC{p) 



(87) 



that is normalized at unity at maximum. The la confidence regions are determined by Sx^ = ~ Xo = 1) while the 
condition 6x'^ = 4 delimited the 2a confidence regions. It is worth stressing that 5x^ = 1 for 1-dimensional likelihoods. 
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Here, Xo is the value of the for the best fit model. Projections of the likelihood function allow to show eventual 
correlations among the model parameters. 

Using the method sketched above, the classes of models which we are going to study can be constrained and selected 
by observations. However, most of the tests recently used to constrain cosmological parameters (such as the SNela 
Hubble diagram and the angular size-redshift) are essentially distance - based methods. The proposal of Dalai et al. 
[127] to use the lookback time to high redshift objects is thus particularly interesting since it rehes on a completely 
different observable. The lookback time is defined as the difference between the present day age of the Universe and 
its age at redshift z and may be computed as : 



where tn = l/J?o = 9-78h~^ Gyr is the Hubble time, and, as above, E{z,p) is the dimcnsionless Hubble parameter, 
where the set of parameters characterizing the cosmological model, {p}, can be taken into account. It is worth 
noticing that, by definition, the lookback time is not sensible to the present day age of the Universe to so that it 
could be possible that a model fits weh the data on the lookback time, but nonetheless it predicts a wrong value for 
to. This latter parameter can be evaluated from Eq.(88) by changing the upper integration limit from z to infinity. 
This shows that it is a different quantity indeed since it depends on the full evolution of the Universe and not only 
on how the Universe evolves from the redshift z to now. That is why this quantity can be explicitly introduced as a 
further constraint. However, it is possible to show from the observations that t^iz) converges to to already at low z 
and then the method can be considered reliable. 

As an example, let us now sketch how to use the lookback time and the age of the Universe to test a given 
cosmological model. To this end, let us consider an object i at redshift z and denote by ti{z) its age defined as the 
difference between the age of the Universe when the object was born, i.e. at the formation redshift zp, and the one 
at z. It is : 



I . _ f°° dz' /■°° _ dz' 



+ z')E{z',p) 

dz' 



{l + z')E{z',p) 



tL{ZF)-tL{z) . (89) 



where, in the last row, we have used the definition (88) of the lookback time. Suppose now we have N objects and 
we have been able to estimate the age ti of the object at redshift Zi for i = 1,2, . . . , N . Using the previous relation, 
we can estimate the lookback time t°^'^{zi) as : 



tt\Zi) = tL{ZF)-t^{z) 

= K'-tiiz)]-K''-tL{zF)] 

= t"^^-ti{z)-Af , (90) 

where t^'^ is the today estimated age of the Universe and a delay factor can be defined as : 



A/ = C - tL{ZF) . (91) 

The delay factor is introduced to take into account our ignorance of the formation redshift zf of the object. Actually, 
what can be measured is the age ti{z) of the object at redshift z. To estimate zf, one should use Eq.(89) assuming a 
background cosmological model. Since our aim is to determine what is the background cosmological model, it is clear 
that we cannot infer Zf from the measured age so that this quantitj^ is completely undetermined. 

It is worth stressing that, in principle. A/ should be different for each object in the sample unless there is a 
theoretical reason to assume the same redshift at the formation of all the objects. If this is indeed the case (as we 
will assume later), then it is computationally convenient to consider A/ rather than Zf as the unknown parameter 
to be determined from the data. Again a likelihood function can be defined as : 



/:,*(p,A/)(xexp[-x?*(p,A/)/2] 



(92) 



with : 
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where Np is the number of parameters of the model, at is the uncertainty on , a,, 



the one on 



i) and the 



superscript theor denotes the predicted values of a given quantity. Note that the delay factor enters the definition of 
xft since it determines t'j^^{zi) from ti{z) in virtue of Eq.(90), but the theoretical lookback time does not depend on 
A/. 

In principle, such a method should work efficiently to discriminate among the various dark energy models. Actually, 
this is not exactly the case due to the paucity of the available data which leads to large uncertainties on the estimated 
parameters. In order to partially alleviate this problem, it is convenient to add further constraints on the models by 
using Gaussian priors'' on the Hubble constant, i.e. redefining the likelihood function as : 



£(p) oc £it(p)exp 
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cxexp[-x'(p)/2] (94) 
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where we have absorbed A/ in the set of parameters p and have defined : 

with h°'"^ the estimated value of h and an its uncertainty. The HST Key project results [109] can be used setting 
(ft, (T/i) = (0.72,0.08). Note that this estimate is independent of the cosmological model since it has been obtained 
from local distance ladder methods. The best fit model parameters p may be obtained by maximizing C{p) which is 
equivalent to minimize the defined in Eq.(95). 

It is worth stressing again that such a function should not be considered as a statistical in the sense that it is 
not forced to be of order 1 for the best fit model to consider a fit as successful. Actually, such an interpretation is 
not possible since the errors on the measured quantities (both ti and to) arc not Gaussian distributed and, moreover, 
there are uncontrolled systematic uncertainties that may also dominate the error budget. 

Nonetheless, a qualitative comparison among different models may be obtained by comparing the values of this 
psoudo oven if this should not be considered as a definitive evidence against a given model. Having more than one 
parameter, one obtains the best fit value of each single parameter pi as the value which maximizes the marginalized 
likelihood for that parameter defined in Eq.(87). After having normalized the marginalized likelihood to 1 at maximum, 
one computes the la and 2a confidence limits (CL) on that parameter by solving Cp- = cxp (—0.5) and Cp, = cxp (—2) 
respectively. In summary, taking into account the above procedures for distance and time measurements, one can 
reasonably constrain a given cosmological model. In any case, the main and obvious issue is to have at disposal 
sufficient and good quality data sets. 



C. Samples of data to constrain models: the case of LSS for lookback time 

In order to apply the method outlined above, we need a set of distant objects whose distances and ages can be 
somehow estimated. As an example for the lookback time method, let us consider the clusters of galaxies which seem 
to be ideal candidates since they can be detected up to high redshift and their redshift, at formation epoch ^ is almost 



* The need for priors to reduce the parameter uncertainties is often advocated for cosmological tests. For instance, in Ref. [128] a strong 
prior on Qm is introduced to constrain the dark energy equation of state. It is likely, that extending the dataset to higher redshifts and 
reducing the uncertainties on the age estimate will allow to avoid resorting to priors. 

^ It is worth stressing that, in literature, the cluster formation redshift is defined as the redshift at which the last episode of star formation 
happened. In this sense, we should modify our definition of A/ by adding a constant term which takes care of how long is the star 
formation process and what is the time elapsed from the beginning of the Universe to the birth of the first cluster of galaxies. For this 
reason, it is still possible to consider the delay factor to be the same for all clusters, but it is not possible to infer zp from the fitted 
value of A/ because we do not know the detail of star formation history. This approach is particular useful since it allows to overcome 
the problem to consider lower limits of the Universe age at z rather than the actual values. 
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Color age 


Scatter age 


z 


N 


Age (Gyr) 


Ref 


2 


N 


Age (Gyr) 


Ref 


0.60 


1 


4.53 


[131] 


0.10 


55 


10.65 


[133] 


0.70 


3 


3.93 


[131] 


0.25 


103 


8.89 


[133] 


0.80 


2 


3.41 


[131] 


1.27 


1 


1.60 


[135] 



Table I: Main properties of the cluster sample used for the analysis. The data in the left part of the Table refers to clusters 
whose age has been estimated from the color of the reddest galaxies (color age), while that of clusters in the right part has 
been obtained by the color scatter (scatter age). For each data point, we give the redshift z, the number A'^ of clusters used, 
the age estimate and the relevant reference. 

the same for all the clusters. Furthermore, it is relatively easy to estimate their age from photometric data only. To 
this end, the color of their component galaxies, in particular the reddest ones, is needed. 

Actually, the stellar populations of the reddest galaxies become redder and redder as they evolve. It is just a matter, 
then, to assume a stellar population synthesis model, and to look at how old the latest episode of star formation should 
be happened in the galaxy past to produce colors as red as the observed ones. This is what is referred to as color 
age. The main limitation of the method rehes in the stellar population synthesis model, and on a few (unknown) 
ingredients (among which tlie nietallirity and the star formation rate). 

The choice of the evolutionary model is a key step in the estimate of the color age and the main source of uncertainty 
[129]. An alternative and more robust route to cluster age is to consider the color scatter (see [130] for an early 
application of this approach). The argument, qualitatively, goes in this way: if galaxies have an extreme similarity 
in their color and nothing is conspiring to make the color scatter surreptitiously small, then the latest episode of star 
formation should happen in the galaxy far past, otherwise the observed color scatter would be larger. 

Quantitatively, the scatter in color should thus be equal to the derivative of color with time multiplied the scatter 
of star formation times. The first quantity may be predicted using population synthesis models and turns out to be 
almost the same for all the evolutionary models thus significantly reducing the systematic uncertainty. We will refer 
to the age estimated by this method as scatter age. The dataset we need to apply the method may be obtained using 
the following procedure. First, for a given redshift Zi, we collect the colors of the reddest galaxies in a cluster at that 
redshift and then use one of the two methods outlined above to determine the color or the scatter age of the cluster. 
If more than one cluster is available at that redshift, we average the results from different clusters in order to reduce 
systematic errors. Having thus obtained ti{zi), we then use Eq.(90) to estimate the value of the lookback time at that 
redshift. 

Actually, what we measure is t°,^^{zi) + dA/ that is the quantity that enters the definition (93) of Xu ^-nd then the 
likelihood function. To estimate the color age, following [131], it is possible to choose, among the various available 
stellar population synthesis models, the Kodama and Arimoto one [132], which, unlike other models, allows a chemical 
evolution neglected elsewhere. This gives us three points on the diagram z vs. t°^'^ obtained by applying the method 
to a set of six clusters at three different redshifts as detailed in Table 1. 

Using a large sample of low redshift SDSS clusters, it is possible to evaluate the scatter age for clusters age at 
z = 0.10 and z = 0.25 [133]. Blakeslee et al. [135] apphed the same method to a single, high redshift {z = 1.27) 
cluster. Collecting the data using both the color age and the scatter age, we end up with a sample of ^ 160 clusters 
at six redshifts (listed in Table 1) which probe the redshift range (0.10,1.27). This nicely overlaps the one probed 
by SNela Hubble diagram so that a comparison among our results and those from SNela is possible. We assume a 
(7=1 Gyr as uncertainty on the cluster age, no matter what is the method used to get that estimate. 

Note that this is a very conservative choice. Actually, if the error on the age were so large, the color - magnitude 
relation for reddest cluster galaxies should have a large scatter that is not observed. We have, however, chosen such 
a large error to take qualitatively into account the systematic uncertainties related to the choice of the evolutionary 
model. 

Finally, we need an estimate of t^^^ to apply the method. Following Rebolo et al. [107], one can choose (ig^%crt) = 
(14.4, 1.4) Gyr as obtained by a combined analysis of the WMAP and VSA data on the CMBR anisotropy spectrum 
and SDSS galaxy clustering. Actually, this estimate is model dependent since Rebolo et al. [107] implicitly assumes 
that the ACDM model is the correct one. However, this value is in perfect agreement with t^'^^ = 12.6^21 Gyr 
determined from globular clusters age [64] and t^^^ > 12.5 ±3.5 Gyr from radioisotopes studies [136]. For this reason, 
one is confident that no systematic error is induced on the adopted method using the Rebolo et al. estimate for tg^* 
even when testing cosmological models other than the ACDM one. 
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D. Dark Energy as a curvature effect 



The methods outHned above allow to constrain dark energy models without considering the nature of dark con- 
stituents. In [98], it is shown that the most popular quintessence (dark energy) models can be reproduced, in principle, 
only considering "curvature effects" i.e. only generalizing the theory of gravity to some f{R) which is not supposed 
to be simply linear in R. From our point of view, this approach seems "economic" and "conservative" and does not 
claim for unknown fundamental ingredients, up to now not detected, in the cosmic fiuid^. As it is clear, from Eq.(20), 
the curvature stress - energy tensor formally plays the role of a further source term in the field equations so that its 
effect is the same as that of an effective fluid of purely geometric origin. Let us rewrite it here for convenience: 

= jj^^ |^5a/3 [f{R) - Rf'iR)] + fiRy^^iga^gp. - ga^g,.)^ ■ (96) 

Our aim is to show that such a quantity provides all the ingredients we need to tackle with the dark side of the 
Universe. In fact, depending on the scales, such a curvature fluid can play, in principle, the role of dark matter and 
dark energy. To be more precise, also the coupling l/f'{R) in front of the matter stress energy tensor, see Eqs.(20), 
plays a fundamental role in the dynamics since it affects, in principle, all the physical processes (e.g. the nucleo- 
synthesis) and the observable (luminous, clustered, baryonic) quantities. This means that the whole problem of the 
dark side of the Universe could be addressed considering a comprehensive theory where the interplay between the 
geometry and the matter has to be reconsidered assuming non-linear contributions and non-minimal couplings in 
curvature invariants. 

From the cosmological point of view, in the standard framework of a spatiafly flat homogenous and isotropic Uni- 
verse, the cosmological dynamics is determined by its energy budget through the Friedmann equations. In particular, 
the cosmic acceleration is achieved when the r.h.s. of the acceleration equation remains positive. Specifically the 
Friedmann equation, in physical units, is 



a 1 

- = - g [Ptot + optot) 



(97) 



The subscript tot denotes the sum of the curvature fluid and the matter contribution to the energy density and 
pressure. From the above relation, the acceleration condition, for a dust dominated model, leads to : 



Pcurv "t" Pm ~t~ '^Pcurv ^ 

SO that a key role is played by the effective quantities : 

Pcurv 

and 



< 



Ptot 



3pc 



{f{R) - Rf'iR)] - SHRf"{R) 



Wr. 



-1 + 



Rf"{R) + R Rf"'{R) - Hf"{R) 



[f{R) - Rf'iR)] /2 - SHRf'iR) 



(98) 



(99) 



(100) 



deduced from Eq.(96). As a flrst simple choice, one may neglect ordinary matter and assume a power -law form 
f{R) = /o-R", with n a real number, which represents a straightforward generalization of Einstein GR in the limit 
n = 1. One can find power -law solutions for a{t) providing a satisfactory fit to the SNela data and a good agreement 
with the estimated age of the Universe in the range 1.366 < n < 1.376 [48, 97]. On the other side, one can develop the 
same analysis in presence of the ordinary matter component, although in such a case, one has to solve numericafly 
the field equations. Then, it is stifl possible to confront the Hubble fiow described by such a model with the Hubble 
diagram of SNela using the above mentioned methods. The data fit turns out to be significant (see Fig. 3) improving 
the value and it fixes the best fit value at n = 3.46 when it is accounted only the baryon contribute flf, « 0.04 
(according with BEN prescriptions). It has to be remarked that considering dark matter does not modify the result 
of the fit, as it is evident from Fig. 3, in some sense positively supporting the assumption of no need for dark matter 



Following the Occam razor prescription: "Entia non sunt multiplicanda praeter necessitatem. " 
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Figure 2: Best fit curve to the SNela Hubble diagram for the power law Lagrangian model. Only data of "Gold" sample of 
SNela have been used. 




Figure 3: The Hubble diagram of 20 radio galaxies together with the "Gold" sample of SNela. in term of the rcdshift as 
suggested in [137]. The best fit curve refers to the i?"- gravity model without dark matter (left), while in the right panel it 
is shown the diflcrcncc between the luminosity distances calculated without dark matter and in presence of this component in 
term of redshift. It is evident that the two behaviors are quite indistinguishable. 

in this model. A part the simplicity of the power law model, the theoretical implications of the best fit values found 
for n are telling us that dynamics related to cosmological constant (whose theoretical shortcomings are well known) 
could be seriously addressed by finding a rehable f{R) gravity model (see also [49]). 

From the evolution of the Hubble parameter in term of redshift, one can even calculate the age of Universe. In Fig. 
4, it is sketched the age of the Universe as a function of the correlation between the deceleration parameter qo and 
the model parameter n. The best fit value n = 3.46 provides tuniv ~ 12.41 Gyr. 

It is worth noticing that considering /(-R) = fo R"" gravity is only the simplest generalization of the Einstein 
theory. In other words, it has to be considered that i?"- gravity represents just a working hypothesis as there is no 
overconfidence that such a model is the correct final gravity theory. In a sense, we want only to suggest that several 
cosmological and astrophysical results can be well interpreted in the realm of a power law extended gravity model. 

As matter of fact, this approach gives no rigidity about the value of the power n, although it would be preferable 
to determine a model capable of working at different scales. Furthermore, we do not expect to be able to reproduce 
the whole cosmological phenomenology by means of a simple power law model, which has been demonstrated to be 
not sufficiently versatile [139-141]. 

For example, we can easily demonstrate that this model fails when it is analyzed with respect to its capability 
of providing the correct evolutionary conditions for the perturbation spectra of matter overdensity. This point is 
typically addressed as one of the most important issues which suggest the need for dark matter. In fact, if one wants 
to discard this component, it is crucial to match the experimental results related to the Large Scale Structure of the 
Universe and the CMBR which show, respectively at late time and at early time, the signature of the initial matter 
spectrum. 

As important remark, we notice that the quantum spectrum of primordial perturbations, which provides the seeds of 
matter perturbations, can be positively recovered in the framework of i?" - gravity. In fact, f{R) oc R? can represent a 
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Figure 4: Contour plot in the plane (go , n) describing the Universe age as induced by i?" - gravity model without dark matter. 
The contours refer to age ranging from 11 Gyr to 16 Gyr from up to down. The dashed curves define the 1 — a region relative 
to the best fit Universe age suggested by the last WMAP release (IS.TStJ'j iy Gyr) in the case of A-CDM model [138]. At the 
best fit n ~ 3.5 for SNcIa, the measured go ~ —0.5 gives a rather short age (about 11.5 Gyr) with respect to the WMAP 
constraint. This is an indication that the /(R) model has to be further improved. 



viable model vifith respect to CMBR data and it is a good candidate for cosmological Inflation (see [144] and references 
therein) . 

In order to develop the matter power spectrum suggested by this model, we resort to the equation for the matter 
contrast obtained in [145] in the case of fourth order gravity (see even [146] for a review on cosmological perturbations 
in f{R) - theories). This equation can be deduced considering the conformal Newtonian gauge for the perturbed metric 
[145] : 

ds"^ = {l + 2ip)dt^ - a'^{l + 2(l))T,^i^^{dx'). (101) 

where ip and </!> are now gravitational perturbation potentials. In GR, it is — —ijj, since there is no anisotropic stress; 
in ETGs, this relation breaks, in general, and the i j components of field equations give new relations between cj) 
and tp. 

In particular, for f{R) gravity, due to the non- vanishing derivatives fn-i-j (with i ^ j), the 4> — ip relation becomes 
scale dependent. Instead of the perturbation equation for the matter contrast 6, we provide here its evolution in term 
of the growth index J-" = dlnS/dlna, which is the directly measured quantity at 2; ~ 0.15 : 



2 



[a) h 



- + -E'{a) 
a a 



jr(a) - - — = , 102 



(the prime, in this case, means the derivative with respect to a, n is the model parameter, being f{R) oc i?"), 
E(a) = H{a)/Ho, R is the dimensionless Ricci scalar, and 

q^_^Jrr^ (103) 

For n = 1 the previous expression gives the ordinary growth index relation for the Cosmological Standard Model. It 
is clear, from Eq.(102), that such a model suggests a scale dependence of the growth index which is contained into 
the corrective term Q so that, when Q ^ 0, this dependence can be reasonably neglected. 

In the most general case, one can resort to the limit aH < k < lO^'^h Mpc^^ , where Eq.(102) is a good approxi- 
mation, and non-linear effects on the matter power spectrum can be neglected. Studying numerically Eq.(102), one 
obtains the growth index evolution in term of the scale factor; for the sake of simplicity, we assume the initial condition 
^{(Hs) = 1 at the last scattering surface as in the case of matter-like domination. The results are summarized in 
Fig. (5) - (6), where we have displayed, in parallel, the growth index evolution in i?" - gravity and in the ACDM model. 
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Figure 5: Scale factor evolution of the growth index : (left) modified gravity, in the case Qm ~ ^bar ~ 0.04, for the SNela 
best fit model with n = 3.46, (right) the same evolution in the case of a ACDM model. In the case of Ji" - gravity it is shown 
also the dependence on the scale k. The three cases k = 0.01, 0.001, 0.0002 have been checked. Only the latter case shows a 
very small deviation from the leading behavior. Clearly, the trend is that the growth law saturates to .7^ = 1 for higher redshifts 
(i.e. a ^ 0.001 to 0.01). This behavior agrees with observations since we know that comparing CMB anisotropics and LSS, we 
need roughly 5 oc a between recombination and z ~ 5 to generate the present LSS from the small fluctuations at recombination 
seen in the CMB. 
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Figure 6: The evolution of the growth index in terms of the scale factor when dark matter is included in the whole energy 
budget. Again, the left plot shows the modified gravity evolution for the SNela best fit model with n = 3.46, while the right 
one refers to ACDM model. 



In the case of Clm = ^bar ^ 0.04, one can observe a strong disagreement between the expected rate of the growth 
index and the behavior induced by power law fourth order gravity models. 

This negative result is evidenced by the predicted value of J^{az = 0.15), which has been observationally estimated 
by the analysis of the correlation function for 220000 galaxies in 2dFGRS dataset sample at the survey effective depth 
z = 0.15. The observational result suggests = 0.58 ±0.11 [147], while our model gives .^(0^ = 0. 15) 0.117 (fc = 
0.01), 0.117 (k = 0.001), 0.122 {k = 0.0002). 

Although this result seems frustrating with respect to the underlying idea to discard the dark matter component 
from the cosmological dynamics, it does not give substantial improvement in the case of i?" - gravity model plus dark 
matter. In fact, as it is possible to observe from Fig. (6), even in this case the growth index prediction is far to be 
in agreement with the ACDM model and again, at the observational scale z = 0.15, there is not enough growth of 
perturbations to match the observed Large Scale Structure. In such a case one obtains: J^{az = o.i5) ~ 0.29 {k = 
0.01), 0.29 {k = 0.001), 0.31 {k = 0.0002), which are quite increased with respect to the previous case but still very 
far from the experimental estimate. 

It is worth noticing that no significant different results are obtained if one varies the power n. Of course in the case 
of n ^ 1, one recovers the standard behavior if a cosmological constant contribution is added. These results seem to 
suggest that an ETG model which considers a simple power law of Ricci scalar, although cosmologically relevant at 
late times, is not viable to describe the evolution of Universe at all scales. 

In other words such a scheme seems too simple to give account of the whole cosmological phenomenology. In fact, 
in [145] a gravity Lagrangian considering an exponential correction to the Ricci scalar, f{R) = R + Aexp{—BR) 
(with A, B two constants), gives a grow factor rate which is in agreement with the observational results at least 
in the dark matter case. To corroborate this point of view, one has to consider that when the choice of f{R) is 
performed starting from observational data (pursuing an inverse approach) as in [98], the reconstructed Lagrangian 
is a non- trivial polynomial in term of the Ricci scalar, as we shall see below. 

A result which directly suggests that the whole cosmological phenomenology can be accounted only by a suitable 
non -trivial function of the Ricci scalar rather than a simple power law function. In this case, cosmological equations. 
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Figure 7: Comparison between predicted and observed values of r = tL{z) + A/ for the best fit ACDM model. Data in Table 
I have been used. 
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Figure 8: Comparison between predicted and observed values of t = tL^z) + Af for the best fit f{R) power-law model as in 
Fig. 3. Data in Table I have been used. Also for this test, it is evident the strict concordance with ACDM model in Fig. 7. 

coming from an f{R) action, can be reduced to a linear third order differential equation for the function f{R(z)), 
where z is the redshift. The Hubble parameter H{z) inferred from the data and the relation between z and R can be 
used to finally work out f{R). 

This scheme provides even another interesting result. Indeed, one may consider the expression for H{z) in a given 
dark energy model as the input for the reconstruction of ,f{R) and thus work out a f{R) theory giving rise to the 
same dynamics as the input model. 

This suggests the intriguing possibility to consider observationally viable dark energy models (such as ACDM 
and quintessence) only as effective parameterizations of the curvature fluid [98, 139]. As matter of fact, the results 
obtained with respect to the study of the matter power spectra in the case of i?" - gravity do not invalidate the whole 
approach, since they can be referred to the too simple form of the model. Similar considerations can be developed for 
cosmological solutions derived in Palatini approach (see [148] for details). 

An important remark is in order at this point. If the power n is not a natural number, i?" models could be not 
analytic for R ^ 0. In this case, the Minkowski space is not a solution and, in general, the post-Minkowskian limit 
of the theory could be bad defined. Actually this is not a true shortcoming if we consider i?"-gravity as a toy model 
for a (still unknown) self-consistent and comprehensive theory working at all scales. 

However, the discussion is not definitely closed since some authors support the point of view that no f{R) theories 
with f = R + aR^\ n ^ 1 can evolve from a matter-dominated epoch a{t) cx t^/'"^ to an accelerated phase [141]. This 
result could be the end of such theories, if the phase space analysis of cosmological solutions is not correctly faced. 

In [142], and recently in [143], it is shown that transient matter-dominated evolutions evolving toward accelerated 
phases are actually possible and the lack of such solutions in [141] depends on an incomplete parameterization of the 
phase space. 

In general, by performing a conformal transformation on a generic f{R) gravity theory, it is possible to achieve, 
in the Einstein frame, dust matter behaviors which are compatible with observational prescriptions. In addition, by 
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exploiting the analogy between the two frames and between modified gravity and scalar-tensor gravity, one can realize 
that physical results, in the two conformally related frames, could be completely different. In other words one can 
pass from a non- phantom phase behavior (Einstein frame) to a phantom regime (Jordan frame) [25]. 

Now, we can suppose to change completely the point of view. In fact, we can rely directly with the Jordan frame 
and we can verify if a dust matter regime is intrinsically compatible with modified gravity. 

As a first example, one can cite the exact solution provided in [48], which has been deduced working only in the 
Jordan Frame (FRW Universe). In particular, one is able to find a power law regime for the scale factor whose rate 
is connected with the power n of the Lagrangian f{R) — foR". 

2n^ — 3n + 1 

In other words, one has a(t) = aot" with a = . Such an exact solution is found out when only 

2 — n 

baryonic matter is considered [153, 154]. It is evident that such a solution allows to obtain an ordinary matter 
behavior (a = 2/3) for given values of the parameter n (i.e. n ~ —0.13, n ~ 1.29). 

Such solutions are nevertheless stable and no transition to acceleration phase then occurs. In general, it is possible 
to show that solutions of the type 

a = aoit - h)^^ , (104) 

where w is the barotropic index of standard perfect fluid, arises as a transient phase, and this phase evolves into an 
accelerated solution representing an attractor for the system [142]. In any case, a single solution exactly matching, 
in sequence, radiation, matter and accelerated phases is unrealistic to be found out in the framework of simple 
/(i?)-power law theories. The discussion can be further extended as fohows. 

Modified gravity can span a wide range of analytic functions of the Ricci scalar where f{R) = foR" only represents 
the simplest choice. In general, one can reverse the perspective and try to derive the form of gravity Lagrangian 
directly from the data or mimicking other cosmological models. 

Such an approach has been developed in [98], and allows to recover modified gravity Lagrangians by the Hubble 
fiow dynamics H{z): in particular, it is possible to show that wide classes of dark energy models worked out in the 
Einstein frame can be consistently reproduced by /(i?)-gravity as quintessence models with exponential potential [99]. 

Clearly the approach works also for the case of coupled quintessence scalar fields. In other words, the dynamics of 
H{z), considered in the Jordan frame, is reconstructed by observational data considered in the Einstein frame then 
assuming one of the two frames as the "physical frame" could be misleading. Here we further develop this approach 
with the aim to show, in general, the viability of f{R) gravity to recover a matter-dominated phase capable of evolving 
in a late accelerating phase. 

From a formal point of view, the reconstruction of the gravity Lagrangian from data is based on the relation which 
expresses the Ricci scalar in terms of the Hubble parameter : 

R=-6(^H + 2H^ + ^^ . (105) 

Now, starting from the above the f{R) field equations (20) one can reconstruct the form of f{R) from the Hubble 
parameter as a function of the redshift z exploiting the relation (105) after this expression has been rewritten in term 
of the redshift itself. 

A key role in this discussion is played by the conservation equation for the curvature and the matter fiuids which, 
in the case of dust matter, (i.e. Pm = 0) gives : 

Pcurv ~t~ 37^(1 -t- Wcurv)Pcurv — y'(^) Pm) 

In particular, one may assume that the matter energy density is conserved : 

Pm = ^MPcTita~^ = '^HIVIm{1 + zf (107) 

with z — 1/a — 1 the redshift (having set a(to) = 1), the matter density parameter (also here, quantities 
labelled with the subscript refers to present day {z = 0) values). Eq.(107) inserted into Eq.(106), allows to write a 
conservation equation for the effective curvature fiuid : 



+ 3H{1 + Wcurv) Pcurv = ^H^flMil + z) 

Rf"{R 



Rf"{R) , . 
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Actually, since the continuity equation and the field equations are not independent [98] , one can reduce to the following 
single equation 



H = - 



1 



2f'{R) 



SH^ftMil + zf + Rf"{R)+ 



+R 



Rfm - Hf'iR)] } , 



(109) 



where all quantities can be expressed in term of redshift by means of the relation 3- = — (1 + z)H—. In particular, 

at az 

for a flat FRW metric, one has : 
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Now, we have all the ingredients to reconstruct the shape of f{R) by data or, in general, by the definition of a suitable 
H{z) viable with respect to observational results. In particular, we can show that a standard matter regime (necessary 
to cluster large scale structure) can arise, in this scheme, before the accelerating phase arises as, for example, in the 
so called quiessence model. 

A quiessence model is based on an ordinary matter fluid plus a cosmological component whose equation of state 
w is constant but can scatter from w = —1. This approach represents the easiest generalization of the cosmological 
constant model, and it has been successfuhy tested against the SNela Hubble diagram and the CMBR anisotropy 
spectrum so that it allows to severely constraint the barotropic index w [100]. 

It is worth noticing that these constraints extend into the region w < —1, therefore models (phantom models) 
violating the weak energy condition are ahowed. From the cosmological dynamics viewpoint, such a model, by 
definition, has to display an evolutionary rate of expansion which moves from the standard matter regime to the 
accelerated behavior in relation to the value of w. In particular, this quantity parameterizes the transition point to 
the accelerated epoch. 

Actually, if it is possible to find out a /(i?)-gravity model compatible with the evolution of the Hubble parameter 
of the quiessence model, this result suggests that modified gravity is compatible with a phase of standard matter 
domination. To be precise, let us consider the Hubble flow deflned by this model, where, as above : 



H{z) = Ho^nMd + + nx{l + ^)3(i+-) (114) 

with fix = (1 — ^m) and w the constant parameter deflning the dark energy barotropic index. This deflnition of the 
Hubble parameter implies: 

R = -2.HI \nM{l + zf + nxil - 3w;)(l + zf^^+'^A . (115) 



The ansatz in Eq.(114) allows to obtain from Eq.(109) a differential relation for f{R{z)) which can be solved numer- 
ically by choosing suitable boundary conditions. In particular we choose : 
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f{z = 0) = f{Ro) = 6H^{1 - Qm) + Ro ■ (118) 

A comment is in order here. We have derived the present day values oidf /dz and d'^f/dz'^ by imposing the consistency 
of the reconstructed f{R) theory with local Solar System tests. One could wonder whether tests on local scales could 
be used to set the boundary conditions for a cosmological problem. It is easy to see that this is indeed meaningful. 

Actually, the isotropy and homogeneity of the Universe ensure that the present day value of a whatever cosmological 
quantity does not depend on where the observer is. As a consequence, hypothetical observers living in the Andromeda 
galaxy and testing gravity in his planetary system should get the same results. As such, the present day values oidf /dz 
and d'^f/dz'^ adopted by these hypothetical observers are the same as those we have used, based on our Solar System 
experiments. Therefore, there is no systematic error induced hy our method of setting the boundary conditions. 

Once one has obtained the numerical solution for ,f{z), inverting again numerically Eq.(115), we may obtain z = z{R) 
and finally get f{R) for several values of w. 

It turns out that f{R) is the same for different models for low values of R and hence of z. This is a consequence of 
the well known degeneracy among different quiessence models at low z that, in the standard analysis, leads to large 
uncertainties on w. This is reflected in the shape of the reconstructed f{R) that is almost w - independent in this 
redshift range. 

An analytic representation of the reconstructed fourth order gravity model, can be obtained considering that the 
following empirical function 



ln(-/) = Zi [ln(-i?)]'^ [l + ln(-i?)]'^ 



(119) 



approximates very well the numerical solution, provided that the parameters {li,l2,h-,li) are suitably chosen for a 
given value of w. For instance, for w = —1 (the cosmological constant) it is : 

{h,l2,h,k) = (2.6693,0.5950,0.0719,-3.0099) . 

At this point, one can wonder if it is possible to improve such a result considering even the radiation, although 
energetically neghgible. Rather than inserting radiation in the (114), a more general approach in this sense is to 
consider the Hubble parameter descending from a unified model like those discussed in [41]. In such a scheme one 
takes into account energy density which scales as : 



p{z) = A 1 + 



l + z 



f3-a r 



1 + 



1+Z 
I + ZJ, 
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(121) 



This model, with the choice (a, /?) = (3, 4), is able to mimic a Universe undergoing first a radiation dominated era (for 
z ^ Zs), then a matter dominated phase (for zi, <^ z Zs) and finally approaching a de Sitter phase with constant 
energy. 

In other words, it works in the way we are asking for. In such a case, the Hubble parameter can be written, in 



natural units, as H = 



- . Pi'-) 



and one can perform the same calculation as in the quiessence case. 



As a final result, it is again possible to find out a suitable /(i?)-gravity model which, for numerical reasons, it is 
preferable to interpolate as f{R)/R: 



^(^) = 1.02 X ^ 



R 



Ro 



-0.04 X ( A)0.3i 



+0.69 X A- 



0.53 



(122) 



where i?o is a normahzation constant. This result once more confutes issues addressing modified gravity as incom- 
patible with structure formation prescriptions. In fact, also in this case, it is straightforward to show that a phase of 
ordinary matter (radiation and dust) domination can be obtained and it is followed by an accelerated phase. 
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Furthermore, several recent studies are pointing out that large scale structure and CMBR anisotropy spectrum arc 
compatible with f{R) gravity as discussed in details in [101, 102] for the metric approach and in [103] for the Palatini 
approach. 

In particular, in [101], it is shown that several classes of f{R) theories can tune the large-angle CMB anisotropy, the 
shape of the linear matter power spectrum, and qualitatively change the correlations between the CMB and galaxy 
surveys. All these phenomena are accessible with current and future data and will soon provide stringent tests for 
such theories at cosmological scales [104]. 



E. The stochastic background of Gravitational Waves "tuned" by f{R) Gravity 

As we have seen, a pragmatic point of view could be to "reconstruct" the suitable theory of gravity starting from 
data. The main issues of this "inverse " approach is matching consistently observations at different scales and taking 
into account wide classes of gravitational theories where "ad hoc" hypotheses are avoided. In principle, as discussed 
in the previous section, the most popular dark energy cosmological models can be achieved by considering f{R) 
gravity without considering unknown ingredients. The main issue to achieve such a goal is to have at disposal suitable 
datasets at every redshift. In particular, this philosophy can be taken into account also for the cosmological stochastic 
background of gravitational waves (GW) which, together with CMBR, would carry, if detected, a huge amount of 
information on the early stages of the Universe evolution [149]. Here, we want to show that cosmological information 
coming from cosmological stochastic background of GWs could constitute a benchmark for cosmological models coming 
from ETGs, in particular for f{R). 

As well known, GWs are perturbations /i^j^ of the metric g^i, which transform as 3-tensors. The GW-equations in 
the transverse-traceless gauge are 

nh{ = . (123) 

Latin indexes run from 1 to 3. Our task is now to derive the analog of Eqs. (123) for a generic f{R). As we have 
seen from conformal transformation, the extra degrees of freedom related to higher order gravity can be recast into a 
scalar field being 



and 



The GW-equation is now 



where 



5m. = e'^'^g^. with e''* = f'{R) . 
R = e-'^'f' {R - 600 - 60;50=*) . 

nh{ = 

□ = e-^-^ (□ + 2^;^V;a) . 



h{=g'^Sgu=e-^f'g'^e^'^Sgu=hi 



(124) 
(125) 
(126) 
(127) 
(128) 



Since no scalar perturbation couples to the tensor part of gravitational waves, we have 



which means that hj is a conformal invariant. As a consequence, the plane- wave amplitudes hj{t) = h{t)el exp{ikmX"^), 
where is the polarization tensor, are the same in both metrics. This fact will assume a key role in the following 
discussion. 

In a FRW background, Eq.(126) becomes 

h+ (SH + 2,pjh + k'^a-'^h = (129) 

being a(t) the scale factor, k the wave number and h the GW amplitude. Solutions arc combinations of Bessel's 
functions. Several mechanisms can be considered for the production of cosmological GWs. In principle, we could seek 
for contributions due to every high-energy process in the early phases of the Universe. 

In the case of inflation, GW-stochastic background is strictly related to dynamics of cosmological model. This 
is the case we are considering here. In particular, one can assume that the main contribution to the stochastic 
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background comes from the amplification of vacuum fluctuations at the transition between the inflationary phase and 
the radiation era. However, we can assume that the GWs generated as zero-point fluctuations during the inflation 
undergo adiabaticaUy damped osciUations (~ l/o) until they reach the Hubble radius H^^ . This is the particle 
horizon for the growth of perturbations. Besides, any previous fluctuation is smoothed away by the inflationary 
expansion. The GWs freeze out for a/k ^ and reenter the radius after the reheating. The reenter in the 
Friedmann era depends on the scale of the GW. After the reenter, GWs can be constrained by the Sachs- Wolfe effect 
on the temperature anisotropy AT/T at the decoupling. More precisely, such fluctuations are degenerated with scalar 
fluctuations, but GWs can, in principle, be measured via B-polarization of the CMB. The measurement is very hard 
to be performed, but many experiments in this direction are presently planned. In any case, AT/T can always be 
used to derive constraints. 

If (f) acts as the inflaton, we have (p <^ H during the inflation. Adopting the conformal time dr] = dt/a, Eq. (129) 
reads 

h" + 2^h' + = (130) 

X 

where x = oe'^. The derivation is now with respect to t]. Inside the radius H~^, we have krj ^ 1. Considering the 
absence of gravitons in the initial vacuum state, we have only negative-frequency modes and then the solution of (130) 
is 

h = k^/^^/2/^^Cexp{-ikr]). (131) 
all 

C is the amplitude parameter. At the flrst horizon crossing {aH = k) the averag ed amplitude = {k/2Txf/^ \h\ of 
the perturbation is 

4 = 2^C- (132) 

When the scale a/k becomes larger than the Hubble radius H^^ , the growing mode of evolution is constant, i.e. 
it is frozen. It can be shown that AT/T < A^, as an upper limit to A^, since other effects can contribute to the 
background anisotropy. From this consideration, it is clear that the only relevant quantity is the initial amplitude 
C in Eq. (131), which is conserved until the reenter. Such an amplitude depends on the fundamental mechanism 
generating perturbations. Inflation gives rise to processes capable of producing perturbations as zero-point energy 
fluctuations. Such a mechanism depends on the gravitational interaction and then (AT/T) could constitute a further 
constraint to select a suitable theory of gravity. Considering a single graviton in the form of a monochromatic wave, 
its zero-point amplitude is derived through the commutation relations: 

[h{t,x), TTh{t,y)]=i5\x-y) (133) 

calculated at a flxed time t, where the amplitude h is the fleld and tt/j is the conjugate momentum operator. Writing 
the Lagrangian for h 

C = ^^r''h,^h.,, (134) 
in the conformal FRW metric g^,^, where the amplitude h is conformally invariant, we obtain 

nn = ^= e^^a'h (135) 
dh ^ ^ 



Eq. (133) becomes 



hit,x),h{y,y) =i^!%— ^ (136) 



and the flelds h and h can be expanded in terms of creation and annihilation operators 

^(*' ^) = / [^1)6-"^^ + h*it)e+"^^] , (137) 

h{t, x) = -^^-^ j d^k [h{t)e-"'^ + h*{t)e+"'^] . (138) 
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Figure 9: Evolution of the GW amplitude for some power-law behaviors of a(t) ^ t" , ~ i*" and f(R) 
time and amplitude strictly depend on the cosmological background giving a "signature" for the model. 
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The commutation relations in conformal time are 



(139) 



From (131) and (132), we obtain C = ^/2iT^He where H and (j) are calculated at the first horizon-crossing and, 
being e^'^ = f'{R), the remarkable relation 

= , ^ , (140) 

holds for a generic f{R) theory at a given k. Clearly the amplitude of GWs produced during inflation depends on the 
theory of gravity which, if different from GR, gives extra degrees of freedom. On the other hand, the Sachs- Wolfe effect 
could constitute a test for gravity at early epochs. This probe could give further constraints on the GW-stochastic 
background, if ETGs are independently probed at other scales. 

In summary, the amplitudes of tensor GWs are conformally invariant and their evolution depends on the cosmological 
background. Such a background is tuned by a conformal scalar field which is not present in the standard GR. Assuming 
that primordial vacuum fluctuations produce stochastic GWs, beside scalar perturbations, kinematical distortions and 
so on, the initial amplitude of these ones is a function of the /(i?)-theory of gravity and then the stochastic background 
can be, in a certain sense "tuned" by the theory. Viceversa, data coming from the Sachs- Wolfe effect could contribute 
to select a suitable f{R) theory which can be consistently matched with other observations. However, further and 
accurate studies are needed in order to test the relation between Sachs- Wolfe effect and .f{R} gravity. This goal could 
be achieved very soon through the forthcoming space (LISA) and ground-based (VIRGO, LIGO) interferometers. 



VI. APPLICATIONS TO GALACTIC DYNAMICS 



The results obtained at cosmological scales motivates further applications of ETGs, in particular of f{R) theories. 
In general, one is wondering whether ETG models, working as dark energy models, can also play a role to explain the 
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dark matter phenomenology at scales of galaxies and clusters of galaxies. 

Several studies have been pursued in this direction [150] but the main goal remains that to seek a unified model 
capable of explain dynamics at every scale without introducing ad hoc components. 



A. Dark matter as a curvature effect: the case of flat rotation curves of LSB galaxies 



It is well known that, in the low energy limit, higher order gravity imphes modified gravitational potentials [62, 151]. 
By considering the case of a pointlike mass m and solving the vacuum field equations for a Schwarzschild - like metric 
[152, 154], one gets as exact solution from a theory f{R) = foR", the modified gravitational potential : 



$(r) = - 



Gm 
2r 



where 



12^2 



1 - V3()/(-' + 12;//' - 83/)- + oO// + 1 



6n2 - 4n + 2 



(141) 



(142) 



which corrects the ordinary Newtonian potential by a power -law term. As we will see, it has to he (3 > and then 
n > 0. In particular, the best fit value will be /3 ~ 0.8 and then n = 3.2. Standard units have been considered here. 
In particular, this correction sets in on scales larger than rc which value depends essentially on the mass of the system 
[154]. This quantity deserves some discussion. As shown in [154], it is derived from the initial conditions of the models 
and it correlates with the core masses of the LSB galaxies which we have taken into account. In some sense, it is 
a sort of further gravitational radius, beside the standard Schwarzschild radius, which rules the central mass of the 
galaxies, and then it is different for different systems. It is interesting to note that, given a generic 2n-order theory of 
gravity, it is possible to find out n characteristic radii [62] and it is intriguing to suspect that they could likely rule 
the structure and the stability of the astrophysical self-gravitating structures [157]. But this is a working hypothesis 
which has to be firmly demonstrated. 

The corrected potential (141) reduces to the standard $ oc 1/r for n = 1 as it can be seen from the relation (142). 
The generalization of Eq.(141) to extended systems is straightforward. Wc simply divide the system in infinitesimal 
mass elements and sum up the potentials generated by each single element. In the continuum limit, we replace the 
sum with an integral over the mass density of the system taking care of eventual symmetries of the mass distribution 
[154]. Once the gravitational potential has been computed, one may evaluate tlie rotation curve vf. {r) and compare 
it with the data. For extended systems, one has typically to resort to numerical techniques, but the main effect may 
be illustrated by the rotation curve for the pointHke case : 



. Gm 
vlir) = 



l+(l-/3)(- 



(143) 



Compared with the Newtonian result v"^ = Gm/r, the corrected rotation curve is modified by the addition of the second 
term in the r.h.s. of Eq.(143). For < /3 < 1, the corrected rotation curve is higher than the Newtonian one. Since 
measurements of spiral galaxies rotation curves signal a circular velocity higher than those which are predicted on the 
basis of the observed luminous mass and the Newtonian potential, the above result suggests the possibility that such a 
modified gravitational potential may fill the gap between theory and observations without the need of additional dark 
matter. It is worth noticing that the corrected rotation curve is asymptotically vanishing as in the Newtonian case, 
while it is usually claimed that observed rotation curves are fiat (i.e., asymptotically constant). Actually, observations 
do not probe Vc up to infinity, but only show that the rotation curve is fiat within the measurement uncertainties up 
to the last measured point. This fact by no way excludes the possibility that Vc goes to zero at infinity. 

In order to observationally check the above result, one can take into account samples of low surface brightness 
(LSB) galaxies with well measured HI + Ha rotation curves extending far beyond the visible edge of the system. 
LSB galaxies are known to be ideal candidates to test dark matter models since, because of their high gas content, 
the rotation curves can be well measured and corrected for possible systematic errors by comparing 21 -cm HI line 
emission with optical Ha and [NIT] data. Moreover, they are supposed to be dark matter dominated so tliat fitting 
their rotation curves without this elusive component could be a strong evidence in favor of any successful alternative 
theory of gravity. The considered sample (Table II) contains 15 LSB galaxies with data on the rotation curve, the 
surface mass density of the gas component and TZ- photometric band, disk photometry extracted from a larger 
sample selected by de Blok & Bosma [155]. We assume the stars are distributed in a thin and circularly symmetric 
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Id 


D 


MO 


Td 


THI 


Mhi 


Type 


UGC 1230 


51 


22.6 


4 


5 


101 


58.0 


Sm 


UGC 1281 


5.5 


22.7 


1 


7 


206 


3.2 


Sdm 


UGC 3137 


18.4 


23.2 


2 





297 


43.6 


Sbc 


UGC 3371 


12.8 


23.3 


3 


1 


188 


12.2 


Im 


UGC 4173 


16.8 


24.3 


4 


5 


178 


21.2 


Im 


UGC 4325 


10.1 


21.6 


1 


6 


142 


7.5 


SAm 


NGC 2366 


3.4 


22.6 


1 


5 


439 


7.3 


IB(s)m 


IC 2233 


10.5 


22.5 


2 


3 


193 


13.6 


SBd 


NGC 3274 


6.7 


20.2 





5 


225 


6.6 


SABd 


NGC 4395 


3.5 


22.2 


2 


3 


527 


9.7 


SAm 


NGC 4455 


6.8 


20.8 





7 


192 


5.4 


SBd 


NGC 5023 


4.8 


20.9 





8 


256 


3.5 


Scd 


DDO 185 


5.1 


23.2 


1 


2 


136 


1.6 


IBm 


DDO 189 


12.6 


22.6 


1 


2 


167 


10.5 


Im 


UGC 10310 


15.6 


22.0 


1 


9 


130 


12.6 


SBm 



Table II: Properties of sample galaxies. Explanation of the columns : name of the galaxy, distance in Mpc; disk central surface 
brightness in the TZ band (corrected for galactic extinction); disk scalelength in kpc; radius at which the gas surface density 
equals 1 M©/pc^ in arcsec; total HI gas mass in 10* M©; Hubble type as reported in the NED database. 



Id 







lo 




C 








xVdof 


^rms 


UGC 1230 





83 


± 





02 


-0 


39 


± 





09 





15 


± 





01 


15.9 


± 0.5 


2.97/8 


0.96 


UGC 1281 





38 


± 





01 


-3 


93 


± 





80 





65 


± 





08 


0.64 


± 0.33 


3.48/21 


1.05 


UGC 3137 





72 


± 





03 


-1 


86 


± 





06 





65 


± 





02 


9.8 


± 0.9 


48.1/26 


1.81 


UGC 3371 





78 


± 





05 


-1 


85 


± 





01 





41 


± 





01 


3.3 


± 0.2 


0.48/15 


1.30 


UGC 4173 





94 


± 





02 


-0 


97 


± 





22 





34 


± 





01 


9.37 


± 0.04 


0.12/10 


0.52 


UGC 4325 





79 


± 





07 


-2 


85 


± 





44 





70 


± 





02 


0.50 


± 0.05 


0.09/13 


1.19 


NGC 2366 





96 


± 





14 


-0 


58 


± 





42 





64 


± 





01 


14.5 


± 0.9 


28.6/25 


1.10 


IC 2233 





42 


± 





01 


-3 


50 


± 





05 





64 


± 





01 


1.29 


± 0.06 


6.1/22 


2.10 


NGC 3274 





71 


± 





03 


-2 


30 


± 





19 





55 


± 





03 


2.3 


± 0.3 


17.6/20 


2.7 


NGC 4395 





13 


± 





02 


-3 


68 


± 





31 





14 


± 





01 


7.6 


± 0.3 


37.7/52 


1.40 


NGC 4455 





87 


± 





05 


-2 


32 


± 





07 





83 


± 





01 


0.42 


± 0.04 


3.3/17 


1.12 


NGC 5023 





81 


± 





02 


-2 


54 


± 





05 





53 


± 





02 


0.91 


± 0.06 


8.9/30 


2.50 


DDO 185 





92 


± 





10 


-2 


75 


± 





35 





90 


± 





03 


0.21 


± 0.07 


5.03/5 


0.81 


DDO 189 





54 


± 





08 


-2 


40 


± 





61 





63 


± 





04 


4.2 


± 0.7 


0.44/8 


1.08 


UGC 10310 





72 


± 





04 


-1 


87 


± 





04 





59 


± 





02 


1.39 


± 0.04 


2.90/13 


1.02 



Table III: Best fit values of the model parameters from minimizing x^(/^j logrc, fg). The values of T*, the /dof are reported 
for the best fit parameters (with dof = N — 3 and N the number of datapoints) and the root of the fit 

residuals. Errors on the fitting parameters and the M/L ratio are estimated through the jacknife method hence do not take 
into account parameter degeneracies [154]. 

disk with surface density S(r) = T*/oe.Tp(— r/r^) where the central surface luminosity Jo and the disk scalelength 
are obtained from fitting to the stellar photometry. The gas surface density has been obtained by interpolating the 
data over the range probed by HI measurements and extrapolated outside this range. 

When fitting to the theoretical rotation curve, there are three quantities to be determined, namely the stellar mass- 
to-light (M/L) ratio, and the theory parameters (/3,rc). It is worth stressing that, while fit results for different 
galaxies should give the same /3, Tc is related to one of the integration constants of the field equations. As such, it is 
not a universal quantity and its value must be set on a galaxy - by - galaxy basis. However, it is expected that galaxies 
having similar properties in terms of mass distribution have similar values of Vc so that the scatter in Tc must reflect 
somewhat the scatter in the terminal circular velocities. In order to match the model with the data, we perform a 
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Figure 10: Best fit theoretical rotation curve superimposed to the data for the LSB galaxy NGC 4455 (left) and NGC 5023 
(right). These two cases are considered to better show the effect of the correction to the Newtonian gravitational potential. 
We report the total rotation curve Vc(r) (solid line), the Newtonian one (short dashed) and the corrected term (long dashed). 



likelihood analysis for each galaxy, using, as fitting parameters /?, logrc (with in kpc) and the gas mass fraction'^ 
fg. As it is evident considering the results from the different fits summarized in Table III, the experimental data are 
successfully fitted by the model. In particular, the best fit range of /3 (/3 = 0.80 ± 0.08), corresponding to i?" gravity 
with 2.3 <n < 5.3 (best fit value n = 3.2), seems well overlaps the above mentioned range of n fitting SNela Hubble 
diagram. 

However, these are only preliminary results which do not completely solve the problem of dark matter in galaxies 
by models coming from ETGs and do not fit the growth of structures. In any case, further evidences on the same 
line of thinking are coming from other samples of galaxies (where also high surface brightness galaxies are considered) 
[156], or from galaxy clusters, where the dark matter range is completely different [157]. 



B. Dark Matter Halos inspired by /(Ji) - gravity 

At this point, it is worth wondering whether a link may be found between f{K) gravity and the standard approach 
based on dark matter haloes since both theories fit equally well the same data. The trait - de - union between these two 
different schemes can be found in the modified gravitational potential which induces a correction to the rotation curve 
in a similar manner as a dark matter halo does. As a matter of fact, it is possible to define an effective dark matter 
halo by imposing that its rotation curve equals the correction term to the Newtonian curve induced by .f{B) gravity. 
Mathematically, one has to split the total rotation curve derived from f{R) gravity as v'^{r) = v1j^{r) + v1^^^^^{r) 
where the second term is the correction. Considering, for simplicity a spherical halo where a thin exponential disk is 
embedded, one can write the total rotation curve as v^ir) = disfcl^') + Z5m('') with v^^iski''') the Newtonian disk 
rotation curve and dm(^) ~ GM£,M{T)/r the dark matter one, MuMir) being its mass distribution. Equating the 
two expressions, we get : 

M.M(,) = M.. (JL) 2^-y(l-q.^Xo(.)-V.(.) _ 
with ri = r/rd, So = TJo, Vdiv) = /o(ry/2)i^o(r//2) x A (77/2)^^1(77/2) and« : 

POO 

Mv,P)= Mv,v',m^-''v'^e-^'dv' (145) 
^0 

with jFo only depending on the geometry of the system and "wir" indicating virial quantities. Eq.(144) defines the 
mass profile of an effective spherically symmetric dark matter halo whose ordinary rotation curve provides the part of 
the corrected disk rotation curve due to the addition of the curvature corrective term to the gravitational potential. It 



This is related to the M/L ratio as T* = [{l- fg)Mg]/{fgLd,) with Mg = IAMhi the gas (HI + He) mass, = T*La and = Stt/ot-^ 
the disk total mass and luminosity. 
* Here Ii and Ki, with I = 1,2 are the Bessel functions of first and second type. 
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Figure 11: Best fit curves superimposed to the data for the total sample of 15 LSB galaxies considered. 



is evident that, from an observational viewpoint, there is no way to discriminate between this dark halo model and a 
f{R) power-law gravity model. Having assumed spherical symmetry for the mass distribution, it is straightforward to 
compute the mass density for the effective dark halo as pDM{r) = {l/A-Kr'^)dMDM /dr. The most interesting features of 
the density profile are its asymptotic behaviors that may be quantified by the logarithmic slope aoM = dhipoM/dhir 
which can be numerically computed as function of i] for fixed values of /? (or n). As expected, aoM depends explicitly 
on /3, while (rc,So,rd) enter indirectly through rjvir- The asymptotic values at the center and at infinity denoted as 
ao and ctoo result particularly interesting. It turns out that ao almost vanishes so that in the innermost regions the 
density is approximately constant. Indeed, ao = is the value corresponding to models having an inner core such 
as the cored isothermal sphere [158] and the Burkert model [159]. Moreover, it is well known that galactic rotation 
curves are typicahy best fitted by cored dark halo models (see, e.g., [160] and references therein). On the other hand, 
the outer asymptotic slope is between —3 and —2, that are values typical of most dark halo models in literature. In 
particular, for /? = 0.80 one finds (a'o,Qoo) = (—0.002, —2.41), which arc quite similar to the value for the Burkert 
model (0, —3). It is worth noticing that the Burkert model has been empirically proposed to provide a good fit to the 
LSB and dwarf galaxies rotation curves. The values of (ao, ciao) we find for our best fit effective dark halo therefore 
suggest a possible theoretical motivation for the Burkert -like models. Due to the construction, the properties of the 
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effective dark matter halo arc closely related to the disk one. As such, we do expect some correlation between the dark 
halo and the disk parameters. To this aim, exploiting the relation between the virial mass and the disk parameters, 
one can obtain a relation for the Newtonian virial velocity Vyir = GMyir/Ryir ■ 

1-/3 1±£ a 

One can numerically check that Eq.(146) may be well approximated as Aid ^ wliich has the same formal structure 
as the baryonic Tully - Fisher (BTF) relation M{, oc V^;^^ with the total (gas + stars) baryonic mass and Vfiat the 
circular velocity on the flat part of the observed rotation curve. In order to test whether the BTF can be explained 
thanks to the effective dark matter halo wc arc proposing, we should look for a relation between Vyir and VjiaL- This is 
not analytically possible since the estimate of Vfiat depends on the peculiarities of the observed rotation curve such as 
how far it extends and the uncertainties on the outermost points. Therefore, for given values of the disk parameters, it 
is possible to simulate theoretical rotation curves for some values of r,; and measure Vfiat finally choosing the fiducial 
value for Tc which gives a value of Vfiat as similar as possible to the measured one. Inserting the relation thus found 
between Vfiat and Vyir into Eq.(146) and averaging over different simulations, one finally gets: 

log Mb = (2.88 ±0.04) log F/iat + (4.14 ±0.09) (147) 

while a direct fit to the observed data gives [161] : 

log Mfc = (2.98 ±0.29) log V/tet + (3.37 ±0.13) . (148) 

The slope of the predicted and observed BTF arc in good agreement leading further support to the f{R) gravity 
model. The zcropoint is markedly different with the predicted one being significantly larger than the observed one, 
but it is worth stressing, however, that both relations fit the data with similar scatter. A discrepancy in the zeropoint 
may be due to the approximate treatment of the effective halo which does not take into account the gas component. 
Neglecting this term, one should increase the effective halo mass and hence Vyir which affects the relation with Vfiat 
leading to a higher than observed zcropoint. Indeed, the larger is Mg/Mj, the more the point deviate from our 
predicted BTF thus confirming our hypothesis. Given this caveat, we may therefore conclude with confidence that 
f{R) gravity offers a theoretical foundation even for the empirically found BTF relation. 

All these results converge toward the picture that data coming from observations at galactic, cxtragalactic and 
cosmological scales could be seriously framed in ETGs without considering huge amounts of dark energy and dark 
matter. 



VII. DISCUSSION AND CONCLUSIONS 



Extended Theories of Gravity can be considered as the natural extension of General Relativity. Also if they are 
not the final theory of gravity at fundamental level (i.e. quantum gravity), they could be a useful approach to 
address several shortcomings of GR. In fact, also at Solar System scales, where GR has been strongly confirmed, some 
conundrums come out as the indications of an apparent, anomalous, long-range acceleration revealed from the data 
analysis of Pioneer 10/11, Galileo, and Ulysses spacecrafts. Such results are difficult to be framed in the standard 
theory of GR and in its low energy limit [60] . 

Furthermore, at galactic scales, huge bulks of dark matter are needed to provide realistic models matching with 
observations. In this case, retaining GR and its low energy limit, implies the introduction of an actually unknown 
ingredient (a huge amount of missing matter) . 

We face a similar situation even at larger scales: clusters of galaxies are gravitationally stable and bound only if 
large amounts of dark matter are supposed in their potential wells. 

Finally, an unknown form of dark energy is required to explain the observed accelerated expansion of cosmic fluid. 
Summarizing, almost 95% of matter-energy content of the Universe is unknown while we can experimentally probe 
only gravity and ordinary (baryonic and radiation) matter. 

Considering another point of view, anomalous acceleration (Solar System), dark matter (galaxies, galaxy clusters 
and clustered structures in general), dark energy (cosmology) could be nothing else but the indications that gravity 
is an interaction depending on the scale and the assumption of a linear Lagrangian density in the Ricci scalar R, as 
the Hilbert-Einstein action, could be too simple for a comprehensive picture at any scale. 

Due to these facts, several motivations suggest to generalize GR by considering gravitational actions where generic 
functions of curvature invariants and scalar flelds are present. This viewpoint is physically motivated by several 
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unification schemes and by field quantization on curved spacetime [8]. Furthermore, it is well known that revisions 
of GR can solve shortcomings at early cosmological epochs (giving rise to suitable inflationary behaviors [13, 15]) 
and explain the today observed accelerated behavior [45, 46]. These results can be achieved in metric and Palatini 
approaches [162-165]. 

In addition, reversing the problem, one can reconstruct the form of the gravity Lagrangian by observational data 
of cosmological relevance through a "back scattering" procedure [98]. 

All these facts suggest that the theory should be more general than the Hnear Hilbert-Einstein one implying that 
extended gravity could be a suitable approach to solve GR shortcomings without introducing mysterious ingredients as 
dark energy and dark matter which seem without explanation at fundamental level. However, changing gravitational 
side could be nothing else but a matter of taste since flnal probes for dark energy and dark matter could come out 
from the forthcoming experiments as LHC. 

Furthermore, in recent papers, some authors have confronted this kind of theories even with the PPN prescriptions 
considering both metric and Palatini approaches. The results seem controversial since in some cases [166] it is 
argued that GR is always valid and there is no room for other theories while other studies [167-169] flnd that 
recent experiments as Cassini and Lunar Laser Ranging allow the possibility that ETGs could be taken into account. 
In particular, it is possible to define generalized PPN-parameters and several ETGs could result compatible with 
experiments in Solar System [54, 92, 168]. 

In principle, any analytic ETGs can be compared with the Hilbert-Einstein Lagrangian provided suitable values of 
the coefficients. This consideration suggests to take into account, as physical theories, functions of the Ricci scalar 
which slightly deviates from GR, i.e. ,f{R) = foR^^~^^^ with e a small parameter which indicates how much the theory 
deviates from GR and then approximate as 

/o|i?|(i+^) ~ /o|i?| (^1 + £ln \R\ + + . . (149) 

Actually, the PPN - Eddington parameters /? and 7 may represent the key parameters to discriminate among 
relativistic theories of gravity. In particular, these quantities should be significatively tested at Solar System scales by 
forthcoming experiments like LATOR [170] while the today available releases arc far, in our opinion, to be conclusive 
in this sense, as a rapid inspection of Table IV suggests. In other words, ETGs cannot be a priori excluded also at 
Solar System scales. 



Mercury Perihelion Shift 


I27-/3- 1| < 3 X 10"^ 


Lunar Laser Ranging 


4/3-7-3 = -(0.7 ± 1) X 10"^ 


Very Long Baseline Interf. 


|7 - Ij = 4 X 10"* 


Cassini Spacecraft 


7 - 1 = (2.1 ±2.3) X 10"^ 



Table IV : A schematic resume of recent experimental constraints on the PPN-parameters. They are the perihelion shift of 
Mercury [171], the Lunar Laser Ranging [173], the upper limit coming from the Very Long Baseline Interferometry [172] and 
the results obtained by the estimate of the Cassini spacecraft delay into the radio waves transmission near the Solar conjunction 
[174]. 

In this paper, we have outlined what one should intend for ETGs in the metric and in the Palatini approach. In 
particular, we have discussed the higher-order and the scalar-tensor theories of gravity showing the relations between 
them and their connection to GR via the conformal transformations. 

In the so called Einstein frame, any ETG can be reduced to the Hilbert-Einstein action plus one or more than 
one scalar field(s). The physical meaning of conformal transformations can be particularly devised in the Palatini 
approach, as discussed in Sec. IV. After, we have discussed some cosmological and astrophysical appHcations of ETGs. 

Although the results outlined are referred to the simplest class of ETGs, power law f{R), they could represent 
an interesting paradigm. Assuming both metric and Palatini approach, it is possible to investigate the viability of 
/(i?) cosmological models. The expansion rate H = d/a may be analytically expressed as a function of the redshift 
z, so that it is possible to contrast the model predictions against the observations. In particular, the SNela Hubble 
diagram, the gas mass fraction in relaxed galaxy clusters, the lookback time to galaxy clusters, and radio galaxies can 
be used to constrain cosmological parameters by distance and time-based methods. 

Also if such models are, up to now, not completely satisfactory to match all the observations, they allow to recover 
accelerated behavior of Hubble fiuid without any unknown form of dark energy. However, the issue of structure 
formation has to be seriously faced in order to understand if such toy models could give rise to a self-consistent 
alternative theory to GR. 
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Furthermore, it is possible to "tune" the stochastic background of GWs and this occurrence could constitute a 
further cosmological test capable of confirming or ruling out ETGs once data from interferometers, like VIRGO, 
LIGO and LISA, wih be available. 

In addition, the modification of the gravitational potential arising as a natural effect in the framework of ETGs 
can represent a fundamental tool to interpret the rotation curves of spiral galaxies. Besides, if one considers the 
model parameters settled by the fit over the observational data on rotation curves, it is possible to construct a 
phenomenological analogous of dark matter halo whose shape is similar to the one of the so called Burkert model. 
Since Burkert's model has been empirically introduced to give account for the dark matter distribution in the case 
of LSB and dwarf galaxies, this result could represent an interesting achievement since it provides a theoretical 
foundation to such a model. 

By investigating the relation between dark halo and the galaxy disk parameters, a relation between Md and Vfiat, 
reproducing the baryonic Tully- Fisher, can be deduced. In fact, exploiting the relation between the virial mass and 
the disk parameters, one obtains a relation for the virial velocity which can be satisfactory approximated as oc V^^^. 
Even such a result seems intriguing since it provides a theoretical interpretation for a phenomenological relation. 

As a matter of fact, although not definitive, these phenomenological issues can represent a viable approach for 
future, more exhaustive investigations of ETGs. In particular, they support the quest for a unified view of the dark 
side of the Universe. In summary, these results seem to motivate a careful search for a fundamental theory of gravity 
capable of explaining the full cosmic dynamics by the only "ingredients" which we can directly and firmly experience, 
namely the background gravity, the baryonic matter, the radiation and also the neutrinos [175]. 



Acknowledgments 

The authors wish to thank G. Allemandi, A. Borowiec, V.F. Cardone, S. Carloni, S. Nojiri, S.D. Odintsov, A. 
Stabile and A. Troisi for discussions, suggestions and results presented in this paper. 



[1] S. Weinberg, Gravitation and Cosmology, Wiley, 1972 New York N.Y. 
[2] A. Guth, Phys. Rev. D 23, 347 (1981). 

[3] I.L. Buchbinder, S.D. Odintsov, and I.L. Shapiro, Effective Action in Quantum Gravity, lOP Publishing (1992) Bristol. 
[4] H. Bondi, Cosmology, Cambridge Univ. Press, Cambridge (1952). 
[5] C. Brans and R.H. Dicke, Phys. Rev. 124, 925 (1961). 

[6] S. Capozziello, R. de Ritis, C. Rubano, and P. Scudellaro, La Rivista del Nuovo Cimento 4 (1996) 1. 
[7] D.W. Sciama, Mon. Not. R. Ast. Soc. 113, 34 (1953). 

[8] N.D. Birrell and P.C.W. Davies, Quantum Fields in Curved Space, Cambridge Univ. Press, Cambridge (1982). 
[9] G. Vilkovisky, Class. Quantum Grav. 9, 895 (1992). 
[10] M. Gasperini and G. Veneziano, Phys. Lett. 277B (1992) 256. 

[11] G. Magnano, M. Ferraris, and M. Prancaviglia, Gen. Relativ. Grav. 19, 465 (1987). 
[12] J. Barrow and A.C. Ottewill, J. Phys. A: Math. Gen. 16, 2757 (1983). 
[13] A. A. Starobinsky, Phys. Lett. 91B (1980) 99. 

[14] ,J.P. Duruisseau and R. Kerner, Gen. Rcl. Grav. 15, 797-807 (1983). 

[15] D. La and P.J. Steinhardt, Phys. Rev. Lett. 62 (1989) 376. 

[16] P. Teyssandier and Ph. Tourrenc, J. Math. Phys. 24, 2793 (1983). 

[17] K. Maeda, Phys. Rev. D 39, 3159 (1989). 

[18] D. Wands, Class. Quantum Grav. 11, 269 (1994). 

[19] S. Capozziello, R. de Ritis, A. A. Marino. Gen. Relativ. Grav. 30, 1247, (1998). 

[20] S. Gottlober, H.-J. Schmidt, and A. A. Starobinsky, Class. Quantum Grav. 7, 893 (1990). 

[21] T.V. Ruzmaikina and A. A. Ruzmaikin, .JETP, 30, 372 (1970). 

[22] L. Amendola, A. Battaglia-Mayer, S. Capozziello, S. Gottlober, V. Miiller, F. Occhionero and H.-J. Schmidt, Class. 

Quantum Grav. 10, L43 (1993). 
[23] A. Battaglia-Mayer and H.-J. Schmidt, Class. Quantum Grav. 10, 2441 (1993). 
[24] H.-J. Schmidt, Class. Quantum Grav. 7, 1023 (1990). 

[25] S. Capozziello, S. Nojiri, and S.D. Odintsov, Phys. Lett. B 634, 93 (2006). 

[26] L. Amendola, S. Capozziello, M. Litterio, F. Occhionero, Phys. Rev. D 45, 417 (1992). 

[27] S. Perlmutter et al., ApJ, 517, 565, 1999; R.A. Knop et al., ApJ, 598, 102, 2003. 

[28] A.G. Riess et al., AJ, 116, 1009, 1998; J.L. Tonry et al., ApJ, 594, 1, 2003. 

[29] P. de Bernardis et al.. Nature, 404, 955, 2000 

[30] R. Stompor et al., ApJ, 561, L7, 2001 

[31] D.N. Spergel et al. ApJS, 148, 175, 2003. 



42 



[32] G. Hinshaw ct al. Ap. J. ApJS, 148, 135 (2003). 
[33] D.N. Spergel et al. astro-ph/0603449 (2006). 
[34] A.G. Ricss ct al., Ap,J, 607, 665, 2004 

[35] v., Sahni, A. Starobinski, Int. J. Mod. Phys. D, 9, 373 (2000). 
[36] T. Padmanabhan, Phys. Rept., 380, 235 (2003). 

[37] E.J. Copeland, M. Sami, S. Tsujikawa, Int. Journ. Mod. Phys. D 15, 1753 (2006). 
[38] A. Kamenshchik, U. Moschella, V. Pasquier, Phys. Lett. B, 511, 265, 2001. 
[39] T. Padmanabhan, Phys. Rev. D, 66, 021301, 2002 

[40] B.A. Bassett, M. Kunz, D. Parkinson, C. Ungarelli, Phys. Rev. D, 68, 043504, 2003 

[41] V.F. Cardone, A. Troisi, S. Capozziello, Phys. Rev. D, 69, 083517, 2004; S. Capozziello, V.F. Cardone, E. Elizalde, S. 

Nojiri, S.D. Odintsov, Phys.Rev.D 73, 043512, 2006 
[42] A. Lue, R. Scoccimarro, G. Starkman, Phys. Rev. D, 69, 044005, 2004 
[43] K. Freese, M. Lewis, Phys. Lett. B, 540, 1, 2002 

[44] G.R. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B, 485, 208, 2000 

[45] S. Capozziello, Int. J. Mod. Phys. D, 11, 483, 2002; S. Nojiri and S.D. Odintsov, Phys. Rev. D, 68, 123512, 2003; 

Carroll, S.M., Duvvuri, V., Trodden, M., Turner, M.S., Phys. Rev. D, 70, 043528, 2004; G. AUemandi, A. Borowiec, M. 

Francavigiia, Phys. Rev. D, 70, 103503, 2004. 
[46] S. Nojiri and S.D. Odintsov, Int. J. Meth. Mod. Phys. 4, 115 (2007). 

[47] S. Capozziello, V.F. Cardone, E. Piedipalumbo, M. Sereno, A. Troisi, Int. J. Mod. Phys. D, 12, 381, 2003 

[48] S. Capozziello, V.F. Cardone, S. Carloni, A. Troisi, Int. J. Mod. Phys. D 12, 1969 (2003). 

[49] A. A. Starobinsky, astro-ph/0706.2041 (2007). 

[50] B. Li and J.D. Barrow, Phys. Rev. D75, 084010 (2007). 

[51] S. Nojiri and S.D. Odintsov, hep-th/06012113 (2006). 

[52] B. Li, J.D. Barrow and D.F. Mota, submitted to Phys. Rev. D, gr-qc/0705.3795. 

[53] G. Chen, B. Rathra, ApJ, 582, 586, 2003; S. Podariu, R.A. Daly, M.P. Mory, B. Rathra, ApJ, 584, 577, 2003. 
[54] CM. Will, Theory and Experiments in Gravitational Physics (1993) Cambridge Univ. Press, Cambridge. 
[55] F.W. Hehl et al.. Rev. Mod. Phys. 48, 393 (1976). 

[56] S. Capozziello, G. Lambiase and C. Stornaiolo, Ann. Phys. (Leipzig), 10, 713 (2001). 

[57] K. Stelle, Gen. Relativ. Grav. 9 (1978) 353. 

[58] R.H. Sanders, Ann. Rev. Astr. Ap. 2 (1990) 1. 

[59] P.D. Mannheim and D. Kazanas, Ap. J. 342 (1989) 635. 

[60] J.D. Anderson, et al., Phys. Rev. D 65 (2002) 082004. 

[61] O. Bertolami, Ch.G. Bohmer, F.S.N. Lobo, (2007) gr - qc/0704.1733 

[62] I. Quant and H.-J. Schmidt, Astron. Nachr. 312 (1991) 97. 

[63] P. Schneider, J. Ehlers, and E.E. Falco, Gravitational Lenses Springer- Verlag (1992) Berlin. 
[64] L.M. Krauss and M. White, Ap. J. 397 (1992) 357. 

[65] Ferraris M, Francavigiia M and Magnano G 1988 Class. Quantum Grav. 5 L95. 

[66] Sokolowski L M 1989 Class. Quantum Grav. 6 2045 

[67] G. Magnano and L.M. SokoLowski, Phys. Rev. D50, 5039 (1994). 

[68] Dicke R H 1962 Phys. Rev. 125 2163 

[69] Damour T and Esposito-Farese G 1992 Class. Quantum Grav. 9 2093. 

[70] S. Capozziello, R. de Ritis, A. A. Marino, Class. Quantum Grav. 14, 3243 (1997). 

[71] V. Faraoni, Cosmology in Scalar-Tensor Gravity, Kluwer Academic, Dordrecht (2004). 

[72] A. Einstein, Sitzung-ber. Preuss. Akad. Wiss., 414 (1925). 

[73] H.A. Buchdahl, J. Phys. A 12 (8) (1979), 1229. 

[74] M. Ferraris, M. Francavigiia, C. Reina, Gen. Relativ. Grav. 14, 243 (1982). 

[75] M. Ferraris, M. Francavigiia and I. Volovich, Class. Quantum Grav. 11, 1505 (1994). 

[76] D.N. VoUick, Phys. Rev. D68, 063510 (2003). 

[77] B.Li and M.C. Chu, Phys. Rev. D74, 104010 (2006). 

[78] B.Li, K,C. Chan, and M.C. Chu, Phys. Rev. D,in press, astro-ph/0610794 (2006). 

[79] G. AUemandi, M. Capone, S. Capozziello, M. Francavigiia, Gen. Relativ. Grav. 38 (2006) 33. 

[80] X. Meng and P. Wang, GRG 36 (8), 1947 (2004); GRG 36 (12), 2673 (2004). 

[81] A. Einstein, Ann. der Physik 49 (1916) 769. 

[82] L.P. Eisenhart, Riemannian Geometry Princeton Univ. Press, Princeton (1955). 
[83] E. Schrodinger, Space-Time Structure, Cambridge Univ. Press, Cambridge (1960). 
[84] T. Levi-Civita, The Absolute Differential Calculus, Blackie and Son, London (1929). 

[85] O. Klein, New Theories in Physics, 77, Intern. Inst, of Intellectual Co-operation, League of Nations (1938). 

[86] E. Cartan, Ann. Ec. Norm. 42 (1925) 17. 

[87] A. Palatini, Rend. Circ. Mat. Palermo 43 (1919), 203. 

[88] V.I. Arnold, Mathematical Methods of Classical Mechanics, Springer- Verlag, Berlin (1978). 

[89] T. Appelquist, A. Chodos, P.G.O. Freund, Modern Kaluza-Klein Theories, Addison- Wesley, Reading (1987). 

[90] M. Kaku, Quantum Field Theory, Oxford Univ. Press, Oxford (1993). 

[91] M.B. Green, J.H. Schwarz, E. Witten, Superstring Theory, Cambridge Univ. Press, Cambridge (1987). 
[92] Damour T., Esposito-Farese G., Class. Quant. Grav. 9, 2093 (1992). 



43 



[93] M. Ferraris, M.Francaviglia, in: Mechanics, Analysis and Geometry: 200 Years after Lagrange; Editor: M. Francaviglia, 
Elsevier Science Publishers B:V:, (1991). 

S. Capozzicllo, R.dc Ritis, C. Rubano, P. Scudcllaro, Int. Journ. Mod. Phys. D 5, 85 (1996). 
S. Capozziello and R. de Ritis, Class. Quantum Grav. 11 (1994) 107. 
V. Sahni and A. Starobinsky, Int. J. Mod. Phys. D9, 373 (2000). 

S. Capozziello, V.F. Cardone, M. Funaro, S. Andreon, Phys. Rev. D 70, (2004) 123501. 
S. Capozziello, V. F. Cardone, A. Troisi, Phys. Rev. D 71, 043503 (2005). 

S. Capozziello, V. F. Cardone, E. Piodipalumbo, C. Rubano, Class. Quant. Grav. 23, 1205 (2006). 

S. Hannestad, E. Mortsell, Phys. Rev. D 66, 0635088 (2002); A. Melchiorri, L. Mersini, C. J. Odman, M. Trodden, Phys. 
Rev. D 68, 043509 (2003); S. Hannestad, E. Mortsell, JCAP 0409, 001 (2004). 
Y. Song, W. Hu, I. Sawicki, Phys. Rev. D 75, 044004 (2007). 
S. Tsujikawa, astro-ph/0705.1032 (2007). 
T. Koivisto, Phys. Rev. D 73, 083517 (2006). 

S. Boughn and R. Crittenden, Nature (London) 427, 45 (2004); P. Fasalba and E. Gaztanaga, MNRAS 350, L37 (2004); 
M.R. Nolta ct al. Ap. J. 608, 10 (2004); P. Zhang, M. Liguori, R. Bean, and S. Dodelson, astro-ph/0704.1932 (2007). 
S. Capozziello Int. Jou. of Geom. Methods in Mod. Phys. 4 (2007) 53. 
M. Tegmark et al., Phys. Rev. D, 69, 103501, 2004 
R. Rebolo et al., astro - ph/0402466, 2004 
R.A. Daly, S.G. Djorgovski, astro - ph/0403664, 2004 
W.L. Freedman et al., Ap,J, 553, 47, 2001. 

V.F. Cardone, S. Capozziello, V. Re, E. Piedipalumbo, A&A, 382, 792, 2002 

R. Saunders et al., MNRAS, 341, 937, 2003; R.W. Schmidt, S.W. Allen, A.C. Fabian, MNRAS, 352, 1413, 2004 
Y. Wang, P. Mukherjee, ApJ, 606, 654, 2004 
Y. Wang, M. Tegmark, Phys. Rev. Lett., 92, 241302, 2004 
W. Hu, N. Sugiyama, ApJ, 471, 542, 1996 

D. Kirkman, D. Tyler, N. Suzuki, J.M. O'Meara, D. Lubin, ApJS, 149, 1, 2003 
D. Eiscnstein et al., ApJ. 633, 560, 2005. 
M.A. Strauss et al., AJ, 124, 1810., 2005. 
SNAP web page : http : / / snap . Ibl . gov 
V. Faraoni, Class. Quantum Grav. 22, 32352 (2005). 
C. Rubano and P. Scudellaro, Gen. Relativ. Grav. 37, 521 (2005). 
S. Sasaki, PASJ, 48, L119, 1996; U. Pen, New Ast., 2, 309, 1997 

S.W. Allen, R.W. Schmidt, A.C. Fabian, MNRAS, 334, Lll, 2002; S.W. Allen, R.W. Schmidt, S. Bridle, MNRAS, 346, 
593, 2003 

S.W. Allen, R.W. Schmidt, H. Ebeling, A.C. Fabian, L. van Speybrock, astro - ph/0405340, 2004 
V. Eke, J.F. Navarro, C.S. Frenk, ApJ, 503, 569, 1998 
W.L. Freedman et al., ApJ, 553, 47, 2001. 
M. Tegmark et al., astro - ph/0310723. 

N. Dalai, K. Abazajian, E. Jenkins, A.V. Manohar, Phys. Rev. Lett., 87, 141302, 2001. 
J.A.S. Lima, J.S. Alcaniz, MNRAS, 317, 893, 2000 
G. Worthey, Ph.D. Thesis California Univ., Santa Cruz, 1992 
R.G. Bower, T. Kodama, A. Terlevich, MNRAS, 299, 1193, 1998 
S. Andreon, C. Lobo, A. lovino, MNRAS, 349, 889, 2004 
T. Kodama, N. Arimoto, 1997, A&A, 320, 41 
S. Andreon, A&A, 409, 37, 2003. 

J. B. Writhe, A. Loeb, and D.G. Barnes ApJ, 634, 715, (2005). 
M. McQuinn et al. ApJ 653, 815 (2006). 
K Sigurdson, A Cooray, Phys. Rev. Lett. 95, 211303 (2005). 
J. P. Blakeslee et al., ApJ, 596, L143, 2003. 
R. Cayrel et al.. Nature, 409, 691, 2001. 
Daly R.A., Djorgovsky S.G., 2004, Astrophys. J., 612, 652 
Spergel D.N., et a., 2006, arXiv: astro-ph/ 0603449 

Capozziello S., Nojiri S., Odintsov S.D., Troisi A., 2006, Phys. Lett. B (2006) 135. 
L. Amendola, D. Polarski, and S. Tsujikawa, Phys. Rev. Lett. 98, 131302 (2007). 
L. Amendola, R. Gannouji, D. Polarski, and S. Tsujikawa, Phys. Rev. D75, 083504 (2007). 
S. Carloni, P. Dunsby, S. Capozziello, A. Troisi, Class. Quant. Grav. 22, 4839 (2005). 
S. Carloni, A. Troisi, P.K.S. Dunsby, gr-qc/0706.0452 (2007). 

Hwang J.C., Noh H., 1996, Phys. Rev. D, 54, 1460; Hwang J.C., Noh H., 2001, Phys. Lett. B, 506, 13. 
Zhang P., 2006, Phys. Rev. D, 73, 123504 

Mukhanov V.F., Feldman H.A., Brandenberger R.H., 1992, Phys. Rept. 215, 203 
Lahav O., et al., 2002, MNRAS, 333, 961. 

S. Capozziello, V.F. Cardone, M. Francaviglia, Gen. Relativ. Grav. 38 (2006) 711. 

S. Capozziello, Ch. Corda, Int. Jou. Mod. Phys. D 15, (2006) 1119; S. Capozziello, Ch. Corda, M. De Laurentis Mod. 
Phys. Lett. A 22, (2007) 1097; S. Capozziello, M. De Laurentis, M. Francaviglia submitted to Int. Jou. Mod. Phys. D 



[94 
[95 
[96 
[97 
[98 
[99 
[100 

101 
102 

103 
104 

105 
106 
107 
108' 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118' 
119 
120 
121 
122 

123 
124 
125 
126 
127 
128^ 
129^ 
130 
131 
132 
133 
134 



135 
136 
137 
138^ 
139 
140 
141 
142 
143 
144 
145 
146' 
147 
148^ 
149^ 



44 



(2007). 

[150] Milgrom M. 1983, Astroph. Journ., 270, 365; Bekenstein J. 2004, Phys. Rev. D, 70, 083509. 
[151] Stcllc K. 1978, Gen. Rcl. Grav., 9, 353 

[152] Capozziello S., Cardone V.F., Carloni S., Troisi A. 2004, Phys. Lett. A, 326, 292 

[153] S. Capozziello, V.F. Cardone, A. Troisi, JCAP 0608 (2006) 1423. 

[154] S. Capozziello, V.F. Cardone, A. Troisi, MNRAS 375 (2007) 1423. 

[155] de Blok W.J.G., Bosma A. 2002, Astron. Astroph., 385, 816 

[156] C. Prigcrio Martins, P. Salucci, astro - ph/07032443. 

[157] S. Capozziello, E. De Filippis, V. Salzano, in preparation. 

[158] Binney J., Tremaine S. 1987, Galactic dynamics, Princeton University Books 

[159] Burkert A. 1995, Astroph. Journ., 447, L25; Burkert A., Silk J. 1997, Astroph. Journ., 488, L55; Borriello A., Salucci P. 

2001, Mont. Not. Roy. Astron. Soc, 323, 285 
[160] Gentile G., Salucci P. 2004, Mon. Not. Roy. Astron. Soc, 351, 953 
[161] McGaugh S.S. 2005, Astroph. Journ., 632, 859. 

[162] Capozziello, S., Carloni, S. and Troisi, A., Rcc. Res. Dcvel. Astronomy. & Astrophys. 1, 625 (2003), astro-ph/0303041. 
[163] AUemandi G., Borowicc A., Francaviglia M., Phys. Rev. D 70, 043524 (2004). 

[164] VoUick D. N., Phys. Rev. D 68, 063510 (2003); Meng X. H., Wang P., Class. Quant. Grav. 20, 4949 (2003); Flanagan 
E. E., Phys. Rev. Lett. 92, 071101 (2004); Flanagan E. E., Class. Quant. Grav. 21, 417 (2004); Meng X. H., Wang P., 
Class. Quant. Grav. 21, 951 (2004); G. M. Kremer and D. S. M. Alves, Phys. Rev. D 70, 023503 (2004). 

[165] T. Multamaki and L Vilja, Phys. Rev. D, 73, 024018, 2006. 

[166] Olmo G.J., Phys. Rev. Lett. 95, 261102 (2005); Olmo G.J., Phys. Rev. D 72, 083505 (2005). 
[167] Capozziello S., Troisi A., Phys. Rev. D 72 (2005). 

[168] S. Capozziello, A. Stabile, A. Troisi Mod. Phys. Lett.A 21 (2006) 2291. 

[169] AUemandi G., Francaviglia M., Ruggiero M., Tartaglia A., Gen. Rel. Grav. 37, 1891 (2005). 
[170] Turyshev S.G., Shao M., K.L. Nordtvedt K.L., pre-print: gr-qc/0601035. 

[171] Shapiro LL, in General Relativity and Gravitation 12, Ashby N., et al, Eds. Cambridge University Press (1993). 

[172] Shapiro S.S., et al, Phys. Rev. Lett. D 92, 121101 (2004). 

[173] Williams J.G., et al, Phys. Rev. D 53, 6730 (1996). 

[174] Bertotti B., less L., Tortora P., Nature 425, 374 (2003). 

[175] A. Capolupo, S. Capozziello, G. Vitiello. Phys. Lett. A 363 (2007) 53. 



